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O \ Abstract 

^ : 

^T* ' We study anomalous dimensions of (super)conformal Wilson operators at weak and strong 

coupling making use of the integrability symmetry on both sides of the gauge/string correspon- 
dence and elucidate the origin of their single- logarithmic behavior for long operators/strings in 
\ the limit of large Lorentz spin. On the gauge theory side, we apply the method of the Baxter 
Q— operator to identify different scaling regimes in the anomalous dimensions in integrable sec- 
tors of (supersymmetric) Yang-Mills theory to one-loop order and determine the values of the 



Lorentz spin at which the logarithmic scaling sets in. We demonstrate that the conventional 
semiclassical approach based on the analysis of the distribution of Bethe roots breaks down in 
this domain. We work out an asymptotic expression for the anomalous dimensions which is valid 
throughout the entire region of variation of the Lorentz spin. On the string theory side, the 
logarithmic scaling occurs when two most distant points of the folded spinning string approach 
the boundary of the AdS space. In terms of the spectral curve for the classical string sigma 
model, the same configuration is described by an elliptic curve with two branching points ap- 
proaching values determined by the square root of the 't Hooft coupling constant. As a result, 
the anomalous dimensions cease to obey the BMN scaling and scale logarithmically with the 
Lorentz spin. 
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1 Introduction 



It is well known that in four-dimensional gauge theories the anomalous dimensions of composite 
Wilson operators carrying a large Lorentz spin scale (at most) logarithmically with the spin. This 
result is just one of the facets of a more general Sudakov phenomenon pQ and it can be traced 
back to the existence of massless particles of spin one in the spectrum - the gauge fields. The 
logarithmic scaling of anomalous dimensions is a universal feature of all gauge theories ranging 
from QCD to the maximally supersymmetric M = 4 Yang-Mills (SYM) theory. In particular, in 
the simplest case of twist-two Wilson operators with large Lorentz spin N ^> 1, the anomalous 
dimension behaves as (in the adjoint representation of the SU(N C ) group) [2] 



where A = ^ M iV c is the 't Hooft coupling constant and r cusp (A) is the so-called cusp anomalous 
dimension r cusp (A) is not universal however and depends on the theory under consideration. 
It has numerous applications in phenomenology of strong interactions and its calculation both 
at weak and strong coupling regimes is one of the long-standing problems in gauge theories. At 
present, the cusp anomaly is known in perturbation theory to the lowest three orders |U |3] and 
there exists a prediction at strong coupling in the Af = 4 SYM theory based on the gauge/string 
correspondence [HI IB] ■ 

The gauge/string correspondence offers a powerful tool to study the dynamics of four- 
dimensional Yang-Mills theories at strong coupling. It establishes a correspondence between Wil- 
son operators in M = 4 SYM theory and certain string excitations on the AdSs x S 5 background 
[TUI Ip]. For operators carrying large quantum numbers (Lorentz spin, isotopic R— charge,...) their 
scaling dimension at strong coupling can be found as an energy of dual (semi) classical string con- 
figurations propagating on the curved space. As was shown in Ref. jB], the operators of twist 
two with large Lorentz spin N are dual to a folded string rotating with the angular momentum 
iV on the AdS3 part of the target space. The resulting expression for the twist-two anomalous 
dimension takes the form with the cusp anomalous dimension at strong coupling given by 



Later, the dual string picture was generalized to Wilson operators of higher twist in the M = 4 
SYM theory carrying both large Lorentz spin and the R— charge For operators built from 
holomorphic scalar fields carrying a unit isotopic charge, the total R— charge equals the twist, L. 
For such operators, in the dual picture the center-of-mass of the string rotates with the angular 
momentum L along a large circle of S 5 . 

The important difference between the operators of twist two, L = 2, and of higher twist, 
L > 3, is that the latter are not uniquely specified by the total Lorentz spin N. More precisely, 
for L > 3 there exist several (superconformal) operators with the same N. These operators mix 
under renormalization and the size of the mixing matrix rapidly increases with L and N. As 
a consequence, the anomalous dimensions of Wilson operators of high twist L > 3 also depend 
on the integers £ = 1,2,... which enumerate eigenvalues of the mixing matrix (=anomalous 
dimensions) for a given Lorentz spin N. For fixed L and large N, possible values of the anomalous 
dimension occupy a band [S|. The lower boundary of the band scales for N — > oo as in (jl.lj) 
while the upper boundary scales as ~ L r cusp (A) In N. This implies that for operators of twist 



7 (A) = 2r cusp (A)lniV + 0(iV ) 



(1.1) 




(1.2) 
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L = 3,4, .. . the minimal anomalous dimension has the same leading asymptotic behavior for 
N — > oo as the twist-two anomalous dimension, Eq. (jl.lj) . This result is rather general and it 
holds true in a generic Yang-Mills theory including QCD and TV = 4 SYM theory. 

In the present paper, we study the properties of the minimal anomalous dimensions on both 
sides of the gauge/string correspondence in the limit of large twist L and Lorentz spin N. On the 
string side, the corresponding single-trace Wilson operators are dual to a folded string spinning 
with large angular momentum N in the AdS3 part of the anti-de Sitter space and boosted along 
a large circle on the sphere with a large angular momentum L jllj . The energy E of this classical 
string configuration defines the leading asymptotics of the anomalous dimension 7(A) = E—L—N 
of the dual Wilson operator in the M = 4 SYM theory in the strong coupling regime in planar 
approximation. The string theory provides a definite prediction for 7(A) as a function of L and 
N. One finds that 7(A) takes different forms in three regimes ITT]: 

• For iV <C L, in the "short" string limit 

m 2 N , . 

7(A) = A— - + ... (1.3) 

• For N ^> L, in the "long" string limit 

for £ str < 1 , 

(1.4) 

for Ur > 1 , 

with the parameter £ str defined as £ str = X\n 2 (N/L)/L 2 . Here the integer m counts the number 
of times the string is folded onto itself. The minimal anomalous dimension corresponds to a 
single-folded string, m = 1. In that case, for £ str 3> 1, the leading asymptotic behavior of 
7(A) does not depend on the twist L and is the same as for L = 2 operators, Eqs. (jl.ljl and 
(II. 2j) . For £ str < 1, the role of the twist L is to create the "BMN window", i.e., a region in 
the parameter space in which the anomalous dimension has an expansion in powers of the BMN 
coupling A' = X/(ttL) 2 [10J. It is believed that the first few terms in the expansion of 7(A) in 
powers of A' should match similar expressions for the anomalous dimensions of Wilson operators 
of twist L and spin N obtained in the M = 4 SYM theory in the weak coupling regime ^2] • 

On the gauge theory side, the calculation of anomalous dimensions of higher twist operators 
with large Lorentz spin turns out to be an extremely nontrivial task in a generic Yang-Mills theory 
even to one-loop order due to a large size of the mixing matrix. The problem can be overcome 
thanks to hidden integrability symmetry of the dilatation operator f3J EJ UHl , which maps the 
one-loop mixing matrix for Wilson operators of twist L belonging to the so-called holomorphic 
SL(2) sector into a Hamiltonian of the Heisenberg magnet of length L and spin s determined by 
the conformal spin of the quantum fields (for a review see Ref. [03 )• This observation allows one 
to calculate the exact eigenspectrum of anomalous dimensions of Wilson operators of arbitrary 
twist and Lorentz spin in integrable sectors of Yang-Mills theory by means of the Quantum 
Inverse Scattering Method |T7j . In the M = 4 SYM theory, the minimal anomalous dimension of 
Wilson operators built from L scalar fields and carrying the Lorentz spin N can be identified to 
one-loop accuracy as the minimal energy in the eigenspectrum of the SL(2) Heisenberg magnet 
of length L and the total spin N + Ls with s = \ [TH]. The gauge/string correspondence suggests 



7(A) 



A m 

Vx 



\n 2 (N/L) + ... 



7T 



mln(N/V\) + ... 
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that the minimal anomalous dimension defined in this way should match the relations ()1.3j) and 
(jl.4j) in the thermodynamic limit L — > oo. 

It follows from (J1.3)) and (|1.4j) that the anomalous dimensions of higher twist operator depend 
on a "hidden" parameter £ str = \ \v?(N/L)/L 2 and their behavior at strong coupling is different 
for £ str < 1 and £ str 1. For £ str 1 the anomalous dimension does not have a perturbative 
expansion in the BMN coupling A' and scales as ~ IniV. On the gauge theory side, previous 
studies jTHj of the Bethe ansatz for the SL(2) spin chain in the thermodynamic limit L ^> 1 
led to the expression for 7(A) which coincides with (jl.3)l in the limit of short strings and with 
the first relation in ()1.4|) in the limit of long strings. They did not reveal however neither any 
trace of the second, logarithmic regime in (|1.4|) . nor appearance of a new parameter similar to 
£st r . This fact is in contradiction with our expectations that the minimal anomalous dimension 
of higher twist operators should scale logarithmically to all loops as iV 00, Eq. p.ljl . The 
goal of the present study is to unravel the logarithmic scaling of the anomalous dimension both 
in the gauge and string theory and to understand the physical meaning of the parameter £ str and 
its counter-part £ on the gauge theory side. 

We shall revisit the calculation of the energy of the SL(2) Heisenberg magnet of spin s — |, 1, | 
using the method of the Baxter Q— operator ^Hj as a main tool and demonstrate that the one- 
loop anomalous dimension in integrable sectors of Yang-Mills theory has the following scaling 
behavior in the thermodynamic limit L — > 00: 



For N < L 



For iV > L 



m 2 N , . 

7(A) = A 47Z^ + ---' (L5) 

A m 2 



\n\N/L) + . . . , forf<l, 



"(A): ( 2 f L (1.0) 

mkiN + ... , for £ > 1 



2vr 2 

depending on the parameter £ = \n(N/L)/L. The minimal anomalous dimension corresponds to 
m = 1. Here s equals the conformal spin of the field entering the Wilson operator, i.e., s = |, 1, | 
for scalar, gaugino and gluon fields, respectively. For scalar operators, Eq. (jl.5|) and the first 
relation in (|1.6J) coincide with similar expressions in ()1.3|) and (jl.4j) . respectively. Notice that the 
anomalous dimension in (jl.6J) does not depend on the spin s for N ^> L which suggests that the 
two regimes in ()1.6|) are universal in all gauge theories. This is indeed the case for £ ^> 1 since 
the coefficient in front of 2 In N coincides with the cusp anomalous dimension at weak coupling, 
Eq. (JTU). 

For N ^> L and £ < 1, the one-loop anomalous dimension exhibits a novel double logarithmic 
behavior (jl.6|) . It was first discovered from the string theory considerations ^T] and was later 
reproduced on the gauge theory side [TH]. A natural question arises whether similar contribu- 
tions arise at higher loops and whether they can be resummed to all loops. The gauge/string 
correspondence suggests that in the M = 4 SYM theory the anomalous dimension in the region 
N > L > 1 admits a BMN-like expansion (with A' = A/(ttL) 2 and £ = \n(N/L)/L) 



7 (A) 



00 , sn 
n=l ^ ' 
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where the coefficient functions c n (£) do not depend on the coupling constant and have the fol- 
lowing asymptotics for £ — > and £ — > oo 



c„(0 = c , n + ci, n £ + C(£ 2 



c n (£) = ©(!/£ 



2n-lN 



with co, n = (— l) n ( — I)"/ 71 " ^° one "l°°P order, Eq. (jl.7|) matches (for m — 1) both relations in 
()1.6j) . In addition, for £ — > oo the coefficient in front of A n in the right-hand side of (jl.7j) scales 
as ~ In N and determines the n— loop correction to the cusp anomalous dimension in the weak 
coupling regime, Eq. In the strong coupling regime, upon the substitution £ = £,ll?/\/~\, 

the perturbative series in p. 7)1 can be resummed to all loops into the following expression (for 
A ^ oo and £ str = A \n 2 (N/L)/L 2 = fixed) 



7(A) 



L E ( A ' ln T c n {0) + ... = L 

n=l ^ ' 



1 + A' In 2 — - 1 

Lj 



;i.9) 



where the ellipsis stands for subleading corrections. For iV — > 00 this relation also reproduces 
the leading asymptotic behavior of the anomalous dimension in the last regime in (|1.4|) and, as 
a consequence, it leads to the expression for the cusp anomalous dimension at strong coupling, 

Eq. dE2D. 

As was already mentioned, the one-loop anomalous dimension (jl.7j) coincides with the energy 
of the SL(2) spin chain of length L in the thermodynamic limit L — > 00. The latter can be found 
within the Bethe Ansatz approach by systematically expanding the energy in powers of 1/L with 
a help of known semiclassical methods j2D|- To leading order of the semiclassical expansion, 
the Bethe roots condense on two symmetric intervals on the real axis [—a, —b] U [b, a], with the 
boundaries a and b being functions of L/N. This leads jTH] to the expression for the energy 
given in (|1.3|) for N <ti L and in the first relation in (jl.4|) for N ^> L. It is believed that 
subleading corrections to the energy are suppressed by powers of 1/L and, therefore, are small. 
We demonstrate that this assumption is only justified for £ = \n(N/L)/L < 1, while for £ ^ 1 
the semiclassical expansion of the energy becomes divergent indicating the change of asymptotic 
behavior of the anomalous dimension, Eq. ([1.4)1 . The reason why the semiclassical expansion 
fails is that the two cuts [—a, —b] and [b, a] collide at the origin, that is b — > for £ — > 00, and 
the Bethe roots have a nonvanishing distribution at the origin. As a consequence, for N ^> L 
the semiclassical corrections to the anomalous dimension run in powers of £. To one-loop order 
they are described in ()1.7j) by the function Cx(£). We argue that the semiclassical series for Ci(£) 
is divergent for £ > 1 and propose an approach which circumvents this difficulty and allows one 
to determine this function for arbitrary £. The resulting expression for the one-loop anomalous 
dimension is valid in the thermodynamic limit throughout the entire interval of N and reproduces 
correct logarithmic behavior for N ^> L, Eq. (11.6)1 . 

On the string theory side, in the dual picture of the folded string spinning in the AdS 3 x S 1 part 
of the target space, the logarithmic behavior of the anomalous dimension at strong coupling is 
associated with the classical string configuration which has two spikes approaching the boundary 
of the AdS space. Thanks to classical integrability of the string equations of motion [21] , the same 
configuration is described by the spectral (elliptic) symmetric curve endowed with a meromorphic 
differential of quasimomentum possessing a double pole at x — ±\/A' and having a prescribed 
asymptotic behavior at the origin and infinity [221 EH 121] • The branching points of the curve, 
±6 str and ±a str , depend on the ratio L/N and the coupling constant A'. We show that for 
N ^> L and £ str < 1 the branching points admit a regular expansion in powers of A' and, as 
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a consequence, the anomalous dimension exhibits the BMN scaling, Eq. (|1.4jl . For £ str ^> 1, 
& str approaches its minimal value x/V so that the inner boundaries of two cuts [— a str , — 6 s tr] and 
[b stT , a str ] coincide with the position of poles of the momentum differential and cannot collide. 
This nonanalyticity manifests itself through the appearance of \f\ prefactor in the logarithmic 
behavior of the anomalous dimension at strong coupling, Eq. (jl.4|) . 

Our consequent presentation is organized as follows. In Sect.|21 we outline a general framework 
for analysis of one- loop anomalous dimensions in the thermodynamic limit. It is based on the 
semiclassical expansion of solutions to the Baxter equation j2Hl I2H1 12Z| • In Sect. El we apply the 
semiclassical approach to determine the minimal anomalous dimension in the thermodynamic 
limit and demonstrate that the semiclassical expansion breaks down for m(iV/L) ^> L due to 
collision of cuts. Then, we present an approach to go consistently beyond the semiclassical 
expansion and use it to describe the minimal anomalous dimension for large Lorentz spin. In 
Sect. El we analyze the asymptotic behavior of the anomalous dimensions at strong coupling 
based on the string sigma model consideration. Section 0] contains concluding remarks. Some 
technical details of our calculations are summarized in the Appendix. 

2 Anomalous dimensions in gauge theory 

Let us start with the calculation of one-loop anomalous dimensions of (super) conformal Wilson 
operators of arbitrary twist L and Lorentz spin N belonging to integrable sectors of (supersym- 
metric) Yang-Mills theories. For quantum fields transforming in the adjoint representation of the 
gauge group, the operators under considerations have the following generic form 

O N>L (0) = Yl c fc ,.. fci tr{^ 1 X(0) J D^ 2 X(0)... J D^X(0)}, (2.1) 

fci+...+fc L =Af 

where X(0) stands for the so-called "good" component of quantum fields of a definite helicity in 
the underlying gauge theory, D + = D^n^ is the covariant derivative projected onto the light-cone, 
n 2 = 0. The expansion coefficients Ck 1 ...k L are fixed from the condition for On,l{0) to have an 
autonomous scale dependence, i.e., Eq. (|2.1|) has to be an eigenstate of the one- loop dilatation 
operator. Integrability allows one to map the one-loop anomalous dimension of the operators 
(12. lj) into energy e of the noncompact SL(2) Heisenberg spin chain of length L and the total spin 

n + Ls dm d US [2EL 

7(A) = ^£ + 0(A 2 ). (2.2) 

Here the (half-) integer s is given by the conformal spin JB| of the quantum field X(0), that is, 
s = 1/2 for scalars, s = 1 for gaugino fields of helicity ±1/2 and s = 3/2 for gauge fields of 
helicity ±1. 

2.1 Exact solution 

Let us first describe the exact solution for the energy e of the SL(2) magnet of length L and 
single-particle spin s in each site. We shall employ the method of the Baxter Q— operator [F3\ 
which proves to be convenient for analyzing various semiclassical limits of e = e(N, L) including 
the limit of the large spin N and length L f2li\ l2*7j . The method relies on the existence of an 
operator Q(u) which acts on the Hilbert space of the SL(2) spin chain and is diagonalized by all 
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eigenstates of the magnet for arbitrary complex parameter u. Discussing the energy spectrum 
it suffices to study the eigenvalues of the Q— operator that we shall denote by Q{u). The same 
function Q(u) determines the wave function of the magnet in the representation of Separated 
Variables 1 - This allows one to analyze Q(u) in the semiclassical limit with the help of the 
WKB machinery well-known from quantum mechanics and, then, determine e(N,L). 
By construction, Q(u) satisfies the second-order finite-difference equation 

(u + is) L Q{u + i) + (u- is) L Q{u -i)= t L {u)Q{u) , (2.3) 

which can be thought of as a Schrodinger equation for a single-particle wave function in the 
representation of Separated Variables jSH]. Here ti{u) is a polynomial in u of degree L with 
coefficients given by conserved charges 

t L (u) = 2u L + q 2 u L - 2 + ... + q L (2.4) 

The lowest integral of motion q 2 is related to the total spin of the SL(2) chain, N + Ls, 

q 2 = -(N + Ls){N + Ls - 1) + Ls{s - 1) (2.5) 

with N = 0, 1, . . .. 

In what follows we shall refer to Eq. ()2.3j) as the Baxter equation. Taken alone, it does not 
fix the function Q(u) and it has to be supplemented by an additional condition that Q(u) has to 
be polynomial in u |31j . Examining the asymptotic behavior of both sides of (J2.3)) for u — > oo, it 
is easy to see that the degree of Q(u) is fixed by the total spin iV and, therefore, up to an overall 
normalization, one can write 

N 

Q{u) = l[(u - X k ) . (2.6) 

k=l 

One substitutes this ansatz into ()2.3|) . takes u = Xk and finds that the roots Ai, . . . , Ajv satisfy 
the Bethe equations 

X k + IS \ L = -pr X k -Xj-i 

Xk-is ~ . \\ X k - Xj + i' 

Solving the Baxter equation ()2.3|) supplemented by ()2.6|) one obtains quantized values of the 
charges 93, . . . , ?l and evaluates the corresponding energy and quasimomentum as jHIj 

e = i{\nQ{is))' - i{\nQ{-is))' , e* e = . (2.8) 

Q{—is) 

Replacing Q{u) by its expression (j2.fi |) one verifies that these relations coincide with those coming 
from the Algebraic Bethe Ansatz jTTj 

N N ■> 

c = Y.^-^ e* = TT^^. (2.9) 

The cyclic symmetry of the single-trace operators (|2.1j) imposes an additional selection rule for 
the eigenstates of the spin magnet, e ld = 1. Equations ()2.8|) and (|2.9j) allow one to calculate the 
energy of the spin chain and, then, obtain the one-loop anomalous dimension of Wilson operators 
(JZH) with a help of (g2J). 

1 For an interpretation of the Q— operator in string theory see Ref. |3U) . 
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2.2 Quasiclassical approach 

Let us examine the Baxter equation ()2.3j) for N + Ls ^> 1. In this limit, the charge q2 takes 
large negative values and one can apply semiclassical techniques j2H| to construct the solution to 
(|2.3jl . To go over to the semiclassical limit, we introduce two scaling parameters 

V = (N + Ls)-\ ( 3 = sLr] = -^_. (2.10) 

By definition, < (3 < 1 with the boundary values corresponding to 

0^0, /3 L -^1. (2.11) 

The parameter rj ^ 1 will play the role of the Planck constant. One rescales the spectral 
parameter as u = x/rj and introduces the eikonal phase (the Hamilton- Jacobi "action" function) 
S(x) as 

f\ \ - 

Q(x/r)) = r]- N exp y~ s ( x )j > S(x) = r]J2\n(x - r]\ k ) . (2.12) 

The energy and the quasimomentum in Eq. ()2.8|) are then given in terms of the function S(x) by 
the following expressions 

£ = i [S'(i(3/L) - S'(-i(3/L)] , e ld = exp \-[S(if3/L) - S{-i(5/L)\ \ . (2.13) 



V 

It also proves convenient to introduce a notation for the "effective potential" 

t(x) = ( V /x) L t L (x/ V ) = 2 + %+ *§ + ... + %, (2.14) 

with q k = q k r] k . 

In the semiclassical approach [23 one assumes that the function S(x) and the integrals 
of motion q k (with k — 3, . . . , L) admit a systematic expansion in powers of r] 

S{x) = S (x)+r ] S 1 (x) + ... , q k = q ( k 0) +vqi 1) + --- ■ (2.15) 

It is tacitly assumed that the expansion of S(x) is convergent and each term is uniformly 
bounded. 2 This leads to the expansion of the effective potential (J2.14j) . r(x) = r (x) +T]Ti(x) + . . . 
with 

1 -~{0) -~{Q) ^(1) ^{1) 

r (x) = 2-4 + % + ... + %, r 1 (x) = % + ... + %, ... (2.16) 

One substitutes (|2.15|) into the Baxter equation ()2.3|) and equates the coefficients in front of 
powers of r\ to get to leading order 

2 cosp(x) = T (x) , p(x) = S'o(x) + — . (2-17) 

In the finite-gap theory, the function p(x) defines the Bloch-Floquet multiplier in an auxiliary 
linear problem for the Baker-Akhiezer function and has the meaning of the (quasi)momentum 

2 As we will show in Sect. [31 this assumption is justified for ln(iV/L) < L and is invalid otherwise. 
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O EH1 EH] ■ For the first subleading term in the semiclassical expansion one finds in a similar 
manner 

$<«) = -?£W) - ^ (nM + *^(*>) • (2-18) 

It is straightforward to derive the subleading terms S' k>2 (x) but we will not need them for our 
purposes. The obtained expressions for the action functions S' (x) and S[(x) depend on yet 
unknown conserved charges g^ 0) and respectively. Quantization conditions for these charges 
follow from the requirement for Q(x/r)), Eq. (j2.12|) . to be a single valued function of x. 

According to (j2.16j) . r (x) is a polynomial of degree L in 1/x. Solving (j2.17j) . one finds that 
the momentum p(x) is, in general, a double-valued function on the complex x— plane with the 
square- root branching points Xj obeying the condition r (xj) = ±2. It is convenient to introduce 
the function y(x) = 2sinp(x) and define a complex curve 



1 -~<0) -~<0) 

y 2 = 4-r 2 (x), ro ( x )=2-- + % + ... + %-. (2.19) 



For arbitrary complex x, except the branching points y(xj) = 0, the relation (|2.19|) defines two 
values for y(x). Then, y(x) being a double- valued function on the complex x— plane, becomes 
a single- valued function on the hyperelliptic Riemann surface defined by the complex curve Tl- 
This surface has a genus L — 2 and is realized by gluing together two copies of the complex 
a;— plane along the cuts running between the branching points X2j-i and x 2 j. 

For the SL(2) magnet the Bethe roots verifying Eq. (J2.7j) take real values only, ImA^ = 0. 
In the semiclassical limit, 77 — > 0, they condense on finite intervals on the real axis where the 
momentum p(x) takes purely imaginary values [21] • In terms of the hyperelliptic curve ()2.19|) . 
this corresponds to y 2 < 0, or 

Tq {x) > 4 , for x e S = [x 2 l-2, X2L-3] U . . . U [x 4 , x 3 ) U [x 2 , x x ] , (2.20) 

where X\ > x 2 > ■ ■ ■ > x 2 l-i and one of the intervals contains the origin. The total number of 
intervals in (J2.2(Jj) equals L — 1 and the end points Xj are just the branching points of the complex 
curve (j2.19|) . r^Xj) = 4. As follows from ()2.19|) . the curve can be parameterized by the set of 
2L — 2 real branching points as 

A 2L ~ 2 

» ! 4n(>-!)' < 2 - 2i > 

The intervals S have the meaning of regions where the classical motion of the system occurs in 
the separated variables 3 . Later on we shall encounter the situation when, say, j th interval shrinks 
into a point, x 2 j = x 2 j-\, so that the motion on this interval is frozen at the classical level. In 
what follows we shall refer to double point. 

The leading term of the semiclassical expansion, So (a;), is determined by the momentum p(x), 
Eq. pTTjl . As follows from its definition pTfjl 



3 In the finite-gap theory P2j, the same intervals have the meaning of forbidden zones in the auxiliary linear 
problem for the Baker- Akhiezer function. 
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The momentum p(x) takes purely imaginary values on the intervals (|2.20J) and its values at 
the end point of the j th interval coincide, p(x2j-i) = p(x2j), with e tp ( X2 ^ = ±1 for r (x2j) = ±2, 
respectively. Continuing p(x) to the complex x— plane one finds that p'(x) is an analytical function 
on the complex plane with cuts running along the intervals (|2.20j) . It defines a meromorphic 
differential on Ti 

t' (x) 

dp = p\x) dx = 0V ' dx . (2.23) 

From (l2~lHJ) and (|2~2^|l one finds r (x) = 2 - 1/x 2 + 0(1/ x 3 ) so that e^°°) = 1 and 

dx 

dp ~ n ) for x — ► oo (2.24) 

x 2 

where '— /+' correspond to the upper/lower sheet of T^. According to the definition (j2.17|) . p(od) 
is defined modulo 2tt. Choosing the normalization condition p(oo) = 0, one finds that at the end 
points of the j th interval in 1)2.20)1 . the momentum takes the values p(x2j) = p(x2j-i) = —Trj for 
X2j > 0. As a consequence, the differential dp satisfies the normalization conditions [33| 132*] 

x 2 j-i r poo i- 

dxp'(x) = - & dp = 0, 2 dxp(x) = - dp = -2izj . (2.25) 

Here in both equations, the integration in the left-hand side goes over the upper sheet of Y^. 
In the right-hand side of the first relation, the differential dp is integrated over the aj— cycle 
encircling the interval [x2j,X2j-i] in the anticlockwise direction. The contour jj in the second 
relation starts on the upper sheet above x = oo crosses the same interval and then goes to infinity 
on the lower sheet (see Fig. HJ. 

The obtained expressions for S' (x) and S[(x), Eqs. (j2.17j) and (j2.18|) . respectively, depend on 
the conserved charges yet to be determined. To work out the quantization conditions for q k 
one examines the first derivative of the eikonal phase (j2.12|) 

- 1 

The discontinuity of S'(x) across the cuts (j2.2Uj) gives the distribution density of rescaled Bethe 
roots f]\k- Assuming that the roots rjXk take finite values for rj — > 0, one finds the behavior of 
S'(x) at infinity on the upper, physical sheet of as 

S'(x) ~ ^ = — S- , for x — oo , (2.27) 

X X 

with rj and j3 defined in ()2.10|) . Replacing S(x) by its semiclassical expansion ()2.15|) and matching 
the coefficients in front of powers of r\ one obtains 4 

So(*)~^, S[(x)~\, ... (2.28) 
X x z 

According to (|2.6|) . the total number of Bethe roots equals N. For a given energy level they 

4 Later on we shall consider solutions to the Baxter equation satisfying Q(u) = Q(—u), or equivalently S(x) = 
S(—x). For such solutions, S' k (x) are odd functions of x and their asymptotics at infinity involves odd powers of 
x only, that is, S' k>1 (x) ~ 1/x 3 . 



9 



K-) 



Figure 1: The definition of the a— cycles and 7— contours on the Riemann surface Ti. The 
dashed lines represent the part of the path on the lower sheet of the surface. 

are distributed on L — 1 intervals ()2.20|) . Denoting by £j the number of Bethe roots on the j th 
interval, one has 

2^jf dxS\x)= V i J = w ^ J - s , (j = l,...,L-l). (2.29) 

The sum of all a— cycles is homologous to zero and, as a consequence, the sum of all a— periods 
is given by the residue of S'(x) at infinity. Together with ()2.27j) this leads to N — l\ + . . . + £l-i- 
Solving the quantization conditions ()2.29|) one can determine quantized values of the conserved 
charges According to ()2.29|) they depend both on the scaling parameter 7] and the set of 
nonnegative integers £i,...,£l-i- Replacing q k in ()2.16|) . ()2.17|) and ()2.18|) by their quantized 
values, one constructs semiclassical expression for S(x) and, then, applies (|2.13|) to determine 
the energy and quasimomentum, 

qk = qk(£i,---,tL-i',T}), e = s(£ 1 ,...,£ L -i;rj) . (2.30) 

Explicit form of these relations for L = 3 can be found in Refs. EH El E]- in particular, 
the quantized values of the energy and conserved charges exhibit remarkable regularity and 
form trajectories. The flow parameter along each trajectory is given by the total spin rj = 
(N + Ls)^ 1 while the integers £1, . . . ,£l-i enumerate the trajectories and encode a nontrivial 
analytic structure in the eigenspectrum. For given N, the total number of trajectories equals the 
number of partitions of N into the sum of L — 1 nonnegative integers £\, . . . , £l-i- This is in a 
perfect agreement with the number of irreducible components entering the tensor product of L 
copies of the SL{2) modules jHJ. 

2.3 Minimal energy trajectory 

In the semiclassical approach described in the previous section, the energy of the spin magnet, 
or equivalently the one-loop anomalous dimension of Wilson operators, is parameterized by the 
complex curve T^, Eq. (J2.19|) . The genus of the curve g = L — 2 is defined by the length of the 
spin chain whereas its moduli are determined by the quantized values of the conserved charges 
qk which depend in their turn on the set of integers £±, . . . , £l-i- Going over to different parts of 
the spectrum amounts to specifying the integers £±, . . . , £l-i- 
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In this paper we are interested in the eigenstates possessing the minimal possible energy for 
a given total spin N. Such eigenstates belong to a particular trajectory to which we shall refer 
as the minimal energy trajectory (see Fig. |U below). To describe it one has first to identify the 
corresponding integers £±, . . . ,£l-i- We remind that the energy ()2.9|) is determined by zeros of 
the function Q(u), or equivalently the Bethe roots. In the thermodynamical limit they condense 
on the intervals 1)2.20)1 . For a given total spin N, the minimal energy is realized when the Bethe 
roots are located on two symmetric cuts most distant from the origin [20], that is, £% = £l-i = 
N/2 and Ej = for j = 2, ... ,L — 2. The corresponding Q(u) should be an even function, 
Q{u) = Q(—u). According to (|2.8j) . (|2.6JI and ()2.3)) . such states automatically satisfy the cyclic 
symmetry condition e td = 1 and possess the following quantum numbers 

N = even , q 2k +i = , (2-31) 

with < k < (L — l)/2. The fact that £j = implies that j th interval does not contain Bethe 
roots and, therefore, it shrinks into a point, x 2 j-i = %2j- From the point of view of separated 
variables, this means that classically all but two collective degrees of freedom are frozen and the 
classical motion is confined to the two intervals with £\ = £l-i = N/2. For the complex curve 
(12.21)1 . this implies that all branching points except four, x\ = —X2L-2 and x 2 = ~x 2 l~3, become 
the double points, x 2 j-i — x 2j, and the curve T L reduces to the elliptic curve (see Eq. (J2.37)) 
below). 

Let us apply the semiclassical approach to obtain the expression for the energy for L ;§> 1 
along the minimal energy trajectory as a function of N. To begin with, we consider the ground 
state of the SL(2) spin chain. It has the total spin N = and is described by a trivial solution 
to the Baxter equation (J2.3)) . Q(u) — 1, or equivalently S(x) =0. From (J2.8j) . the corresponding 
energy is e = and the integrals of motion can be read off from ()2.3)) upon substituting Q(u) = 1. 
To leading order in rj = l/(sL) the transfer matrix and momentum, Eq. (J2.17)) . look like 

t (x) = 2cos(l/a;), p{x) = l/x. (2.32) 

Matching these expressions into (J2.19)) and ()2.21)) one finds that for L — > 00 all branching points 
of the spectral curve ()2.21)) are double points, x 2 j-i = x 2 j = 1/(717). 

Let us now consider the minimal energy eigenstates with N ^> 1 and L ^> 1 and distinguish 
two limiting cases: (i) N/L = fixed and (ii) N/L ^> 1, corresponding to < (3 < 1 and 
/3 — > 0, respectively, Eq. (J2.10)) . We recall that the function Q(u) has exactly N roots, Eq. (J2.6j) . 
Assuming that the Bethe roots scale as ~ I/77, one finds from (J2.9)) that for N/L = fixed 
the energy should behave as e ~ N/L 2 ~ 1/L. Indeed, one can obtain the same scaling of e by 
naively expanding the energy ()2.13j) in powers of 1/L 

e = <?"(0) + 0(f3 3 /L*) = K(0) + vS'KO) + Otf)] . (2.33) 

Here in the second relation we replaced S(x) by its semiclassical expansion (|2.15|) . In a similar 
manner, one obtains for the quasimomentum ()2.13|) 

6 = 2s [S' Q (0) + T}S[(0) + 0( V 2 )] = , (2.34) 

where we took into account that Sk(x) are even functions of x for the minimal energy states 

So(0) = 0, S[(0) = 0, .... (2.35) 
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Then, one finds from 1)2.17)1 the asymptotic behavior of the momentum around the origin on the 
upper ( "physical" ) sheet of the Riemann surface ()2.19j) as 

p(x) = ^ + O(x) , (2.36) 

so that the differential dp has a double pole above x = 0. It is important to stress that the 
relations (|2.33|) and (|2.36|) were obtained under the assumption that S'(x) is regular at the 
origin. For the minimal energy eigenstates, all but two intervals in ()2.20|) shrink into the double 
points and the Bethe roots condense on two symmetric cuts on the real axis that we shall denote 
as [—a, —b] and [b, a]. Then, the spectral curve ()2.21|) reduces to 



2 ~2 

y =y 



2 . 12 



y 2 = (x 2 - a 2 )(x 2 - b 2 ) , (2.37) 



where we took into account that for the minimal energy states the branching points appear 
in pairs of opposite sign. The parameters a, b and x<ij are determined by the condition for 
S' (x), Eq. (J2.17j) . to be an analytical function on the complex a;— plane with two symmetric 
cuts on the real axis and prescribed asymptotic behavior at x = and x — > oo, Eqs. ()2.35j) 
and (j2.28|) . respectively [3SJI32]- Equation ()2.33|) implies that the energy has the BMN scaling 
in the thermodynamic limit, e ~ 1/L. This relation is in an apparent contradiction with the 
well-known fact that for N ^> L the anomalous dimension should scale logarithmically e ~ In N 
with the prefactor being L independent. To understand the reason for this discrepancy, one 
notices that the relation ()2.33j) is valid provided that S'(x) is analytical in the vicinity of x = 0, 
or equivalently, there is no accumulation of the Bethe roots around the origin. For N, L — > oo 
with N/L = fixed this is indeed the case but, as we will argue in the next section, the situation 
drastically changes for N ^> L. 

As a hint, let us consider the eigenstates with L = fixed, N —>■ oo and the charges satisfying 
(I2.31|) . In this case, the transfer matrix t (x), Eq. (|2.16j) . is an even polynomial of degree L 
which scales for x — > as 

^<0) 

r (x) ~ \ (2.38) 
x^ 

provided that q£ ^ 0. If g^ 0) vanishes, the asymptotics of tq(x) is governed by the first non- 
vanishing charge q^l ^ with 2k < L. Substituting (|2.38j) into (j2.23j) one finds the leading 
asymptotic behavior of the momentum for x — > on the upper sheet of ()2.19|) as 

fl T 

dp ~ iL— , (2.39) 

x 

yielding p(x) ~ iL In x. According to ()2.10|) . (3 ~ sL/N — > for N — > oo and one finds from 
f!2.17|) the asymptotics of S' Q (x) above i = 0on the physical sheet of (|2.19|) as 

S' Q (x) ~ iLhix . (2.40) 

Comparing this relation with (J2.35j) . one concludes that S' (x) is no longer analytical at the origin 
due to accumulation of rescaled Bethe roots at x — 0. Finally, one substitutes (|2.4(J|) into (|2.13j) 
and obtains the energy 

e ~ 2L \n(L/f3) ~ 2L In iV . (2.41) 
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Notice that the coefficient in front of In N in the right-hand side of 1)2.41)1 is determined by the 
leading asymptotic behavior of the transfer matrix (J2.38)) for x —* 0. For q^ 1 ^ this coefficient 
takes the maximal possible value. For q£ = . . . = = 0> the transfer matrix (|2.1fij) scales as 

To 0*0 ~ Q2m/ x2m (recall that the charges with odd indices vanish, Eq. (|2.31j) ) leading to 

e~4mlniV, ^ +2 = ... = ^) = 0. (2.42) 

The values of the remaining charges qf\ . . . , w ^ be determined below (see Eq. (12.48)1 ). 
Throughout this paper we are interested in the eigenstates with the minimal energy for given N. 
Obviously, they correspond to m — 1, that is, 



e 



(-=i) ~ 4 in N , qf ] = . . . = q$> = . (2.43) 



Here the superscript '(o)' refers to the leading order approximation of the semiclassical expansion 
and q^l = does not necessarily imply that %k = but rather g 2 fc = 0(v)i or equivalently 
<? 2 fc < (-<? 2 ) fc . 

The transfer matrix tq(x), Eq. ()2.1(i)l . corresponding to the minimal energy state ()2.43j) is 
given by 



4 m - 1 \x) = 2-- (2.44) 



1 



and the spectral curve ()2.19|) looks like 



If": y 2 ={ x2 -^)^- < 2 ' 45 > 

It is easy to see that r^ -1 '* coincides with the complex curve for the spin chain of length 2, that 
is Tl=2, Eq. ()2.19j) . Indeed, choosing qfm+i = • • • = Q^l = in ()2.19|1 one effectively descends 
from an infinitely long spin chain (for L — > oo) to the one with finite length 2m. Comparing 
()2.45|) with ()2.21|) we conclude that for m = 1 all but two branching points condense at the 
origin, b = X2j = 0, and the two remaining (resolved) branching points are located at ±1/2. For 
m > 2, one can show [HB] that the complex curve corresponding to ()2.42|) is given by 



with l/x2j = 2m cos (^-)- Comparison with (|2.37p yields a = l/(2m) and 6 = 0, that is, the two 

cuts [—a, —b] and [b, a] collide at the origin. Notice that r^ 11 " 2 "* has exactly 2m double points 
satisfying a;|- > l/(2m) 2 while the remaining double points in ()2.37|) belonging to the interval 
\—b, b] condensed at x = as b — > 0. 

The complex curves ()2.45|) and ()2.46|) have genus zero and, as a consequence, the momentum 
p'(x) is expressed in terms of elementary functions. Indeed, one substitutes ()2.45|) and ()2.46|) 
into ()2.23|) and obtains the following expression for the momentum for arbitrary m 

dx r ( 1 \ 

dp = — , p(x) = dp = 2marcsin I . (2.47) 

1 _ V 2mx / 



x 2 — 



{2m) 2 
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Figure 2: Symmetric two-cut configuration (a) resulting in the BMN scaling of the anomalous 
dimension in gauge theory for £ < 1. For £ 3> 1 the two cuts collide at the origin (b) yielding the 
logarithmic scaling. The same configuration in string sigma model (c) for £ str ^> 1 - the minimal 
value for the inner end of the cut is given by the BMN coupling V^V which prevents the cuts to 
collide. 



Together with (|2.17|) and ([2.16)1 this leads to the following expressions for the integrals of motion 

_ 2(-iyw- 2 "r(m + n) 

q2n ~ r(2n + l)r(m-n + l)" 1 J 

We conclude from ()2.47|) that for m — 1, for the eigenstate with the minimal energy, in the limit 
N —>■ oo and L = fixed, the momentum p'(x) is an analytical function on the complex plane 
with the square root cut [— |, |] and a simple pole at the origin. This should be compared with 
the analytical properties of the momentum in the region N/L = fixed in which case p'(x) has a 
double pole at x — 0, Eq. (|2.36j) . 

We demonstrated in this section that the analytical properties of the momentum p(x) and the 
action function S'(x) are quite different in the two limits mentioned above. For < (3 < 1, the 
Bethe roots condense on two symmetric intervals [—a, —6] and [b, a] with the end-points a and b 
depending on (3. For (3 — > 0, one has b — > and a — > | so that the two cuts collide and form a 
single cut [— |, |] (see Fig. El (a) and (b)). This leads to different scaling behavior of the energy 
in these two limits, Eqs. ()2.33|) and ()2.43|) . respectively. The question remains what happens 
in the intermediate region of the parameter (3 = Ls/(N + Ls), Eq. (j2.10|) . and how important 
the corrections ~ r] to the energy ()2.33|) are. We shall demonstrate in the next section that the 
semiclassical expansion of the energy (J2.33)) becomes divergent for (3 — > due to the collision of 
cuts and work out an asymptotic expression for the energy valid for small (3. 



3 Symmetric two-cut solution in thermodynamic limit 

In this section, we construct the minimal energy trajectory in the thermodynamic limit L 3> 1 
for arbitrary values of the spin N. To begin with, we work out the semiclassical expansion of 
solutions Q(u) to the Baxter equation ()2.3|) for N/L = fixed and, then, extend the analysis to 
the region N 3> L. 



3.1 Semiclassical expansion at leading order 

To leading order of the semiclassical expansion (j2.15|) and (j2.12|) . the function S' (x) is determined 
by the momentum p(x), Eq. f)2 . 1 7j) . The latter is fixed by the condition that dp(x) should be 
a meromorphic differential on the complex curve 1)2.37)1 with prescribed asymptotic behavior at 
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infinity, Eq. (|2.24J) . and at the origin, Eq. ()2.36|) . on the upper sheet of Tl 

(3.1) 



dp ~ 9 , dp ~ dx 



X 2 



4 + O(x ) 



with /3 = Ls/(N + Ls), Eq. (|2.1()j) . Being combined together, these conditions lead to the 
following expression [21] 

dp = x dx . (3.2) 
a/ (a; 2 — a 2 )(x 2 — b 2 ) 

It defines dp as a meromorphic differential on the elliptic curve y 2 = (x 2 — a 2 )(x 2 — b 2 ), Eq. (|2.37|) . 
which has a pair of double poles above x = on both sheets of the corresponding Riemann surface. 
The values of the parameters a and b are fixed by the normalization conditions ()2.25|) 

POO 

dp = , / dp = —nm . (3-3) 

J a 

Here the integer m defines the position of the interval [b, a] inside (|2.2(Jj) . We remind that in 
the limit we are interested in, the complex curve Tl gets reduced to the genus one curve ()2.37|) . 
Namely, all cuts in ()2.2U|) except two, [—a, —b] and [b, a], shrink into points. The integer m — 1 
counts how many collapsed cuts are situated to the right from the interval [b, a] on the real axis. 
The energy depends on m and, as we will see in a moment, it takes the minimal value for m = 1, 
that is, when all shrunken cuts are located inside the interval [—b,b\. 

Solving (|3.3j) one finds the positions of the cuts as a function of the parameter ft, Eq. (I2.10J) 

1 p 1 



" 2toE(t)' " 2mK(r)' ^ ^E(r) ' (3 ' 4) 

where IK(r) and E(r) are elliptic integrals of the first and second kind, respectively, and the 
modular parameter is defined as (for b < a) 

b 2 

r = l- — . (3.5) 

cr 

From ()3.2j) and ()2.17j) one obtains the leading term in the semiclassical expansion as 

n I i gab n 

S%(x) = p\x) + 4 = x2 = + 4 > ( 3 - 6 ) 

0V KJ x 2 ^(x 2 -a 2 )(x 2 -6 2 ) x 2 ' 

with the parameters a and b given by ()3.4|) . One verifies that Sq(x) has a regular expansion 
around x = which is in agreement with our expectations that the Bethe roots condense on the 
intervals [—a, —b] U [b, a] and there is no accumulation of roots at the origin. The energy ()2.33|) 
is determined by the leading term in the expansion of this expression ^H] 



2/3 

£ = ~T 



1 /3 ( 1 

+ 



ab 2 V a 2 b 2 



^K(r) [(2 - r)K(r) - 2E(r)] . (3.7) 



The dependence of e on the total spin iV and conformal spin s enters into this expression through 
the parametric dependence of the modular parameter r, Eq. ()3.4|) . on the scaling parameter 
(3 = Ls/(N + Ls), Eq. (|2~TU|) . 
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We observe that the minimal energy corresponds to m = 1, that is, when all double points 
X2j, Eq. ()2.37|) . are located inside the interval [—b,b]. It is straightforward to find their position 
on the real axis with the help of (|2.32|) and (|3.2j) . Since the total number of double points is 
infinite for L —>■ oo and they occupy a compact interval, x 2 j should be a smooth function of j. 
Differentiating the second relation in ()2.32|) with respect to j, one finds 



dx 2j * = 2 f, , 1( .2 

dj P'(x 2j ) - 23 (3 



[1 + 0{x%)] , (3.* 



where we substituted the momentum by its value at the origin (|3.1|) . It follows from ()3.8|) that 
for large j the double points l/x 2 j are equidistantly distributed on the real axis (see Fig. OJ) 

— = lj + 0(1/ j) ■ (3-9) 

x 2j p 

We remind that the double points x 2 j verify the relation t (x 2 j) = ±2. It is instructive to 
compare (j3.9|) with a similar relation l/x 2 j = nj for the momentum ()2.32|) corresponding to the 
state with N = 0, or equivalently (3 = 1. Equation ()3.9j) suggests that close to the origin and 
away from the cuts, x 2 <C b 2 , the transfer matrix for the two-cut solution looks like r (x) = 
2cosp(x) ~ 2cos(/9/x). Matching this relation into the general expression for t (x), Eq. ()2.14|) . 
we find that the integrals of motion scale as 

V 2n ~ 2 j2^\ p2n > ( 3 - 10 ) 

where q 2n = q 2n /(N + Ls) 2n . For (3 — 1, or equivalently N = 0, this relation is exact for all q 2 k- 
For (3 — > 0, or equivalently iV — - > oo, it cannot be exact since, by definition, q 2 = — 1 + 0{rj). It 
instead indicates that higher charges take anomalously small values q 2 k = 0(rj 2k ) with rj ~ 1/N. 
This is in agreement with our expectations that for (3 = the transfer matrix reduces to ()2.44|) . 

Let us examine the expression for the energy ()3.7|) in two limiting cases /3 — > 1 and (3 — > 0, or 
equivalently N <^ L and N ^ L, respectively, Eq. (|2.11j) . According to (|3.4jl . the corresponding 
values of the modular parameter are r — >• and r — > 1. For r — > 0, one finds from (|3.4jl 

a = — (l + j) +-.., 6 = — (l- /3 = l-^ + ..., (3.11) 
nm V 4/ 7rm V 4/ 16 

where the ellipses denote sub leading terms. Since a — b ~ r — > 0, the two cuts shrink into points 
as /3 — > 1. One expands ()3.7|) in powers of r and obtains the energy as 

l»jv m 2 , o iV 27r 2 m 2 

For r — ► 1, one has from (|3.4jl 

a=— + ..., 6=^-^ + ..., /3 = v^T^ln— = + (3.13) 
2m 2m ^/l — r 

According to ()2.10j) . the limit (3 — > corresponds to /3 = Ls/N with N ^> L. Then, making use 
of ()3.7)1 one finds for the energy 

L«iv 4m 2 2 VI - r Am 2 2 N 

e = —In — 4 — + ... = -^-1x1 - + .... (3.14) 

This relation matches the expression for the energy of long spinning string folded m times [TJ 
and agrees with the thermodynamic Bethe ansatz analysis 
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Figure 3: The transfer matrix cosp(x) = r (x)/2 as a function of 1/x for the symmetric two-cut 
solution for j3 = 3/4 and m = 1. The cuts [—a, —6] U [6, a] with a = 0.45 and b = 0.15 correspond 
to cos 2 p(x) > 1. The double points x 2 j are denoted by crosses, cosp(x 2 j) = ±1 and l/xfj > 1/b 2 . 
The "large" and "small" roots of the transfer matrix, cosp(<5 n ) = 0, are shown by full and light 
blobs, respectively. 

3.2 Breakdown of semiclassical expansion 

For m = 1 the expressions ()3.7|) . (|3.12|) and ()3.14|) describe the dependence of the energy on 
the total spin N + Ls along the minimal energy trajectory at leading order of the semiclassical 
expansion for L ^> 1 and arbitrary values of N including the regions iV <C L and N ^> L. From 
fl2.33|) one expects that the contribution of the subleading terms to the energy is suppressed by 
powers of rj = 1/(L + Ns). Let us compare these predictions with the exact expressions for the 
minimal energy obtained from numerical solution to the Baxter equation, Eqs. ()2.3|) and ()2.8|) . 
For L = 10, the comparison is shown in Fig. |1] on the right panel. We observe that the curve 
determined by ()3.7|) agrees with the exact numbers for N < L and deviates significantly from 
them for N > L. Most importantly, forN^>L the exact energy scales as ~ 41niV while the 
semiclassical formula (J3.14j) yields a different asymptotic behavior. 

To understand the reason for this discrepancy, let us revisit the calculation of the energy 
in the previous section and pay special attention to the position of the cuts. By construction, 
the cuts run along two symmetric intervals [—a, —b] U [b, a}. According to ()3.13|) . for r — > 1 the 
parameter a approaches the value a = l/(2m) while the parameter b vanishes indicating that the 
two cuts collide at the origin. Then, one examines the expression for the momentum (|3.2j) and 
finds that for b = the differential dp reduces to (|2.47|) . We argued in Sect. 12.31 that Eq. (|2.47j) 
leads to a new, logarithmic scaling of the energy ()2.43|) . This suggests that the problem with 
recovering this regime within the semiclassical approach does not lie in the construction of the 
momentum p(x) but rather in the expression for the energy ()2.33|) . 

The exact formula for the energy (J2.8|) involves S'(x±) evaluated at x± = ±i/3/L. The relation 
(12.33)) was derived under the assumption that S'(x) is analytical at the vicinity of these points. 
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For the two-cut solution constructed in the previous section this assumption is justified provided 
that the inner boundary of the cut b stays finite in the limit L — > oo. For x± — ► with b = fixed, 
the semiclassical analysis is applicable and one arrives at the asymptotic behavior (|3.12j) and 
(j3.14j) depending on the ratio N/L. The situation is more complex when b vanishes in the scaling 
limit. For b — > with x± = fixed, the two cuts collide and the very assumption that S'(x) is 
analytical at the origin does not hold true anymore. In this case, one expects that the energy 
will develop a new asymptotic behavior and the crossover will occur for b ~ |x±| (see Fig. El (a) 
and 0>)). 

Defining a new parameter £ = |x±|/£> = (3/(bL) one anticipates that the semiclassical approach 
is applicable for £ < 1. Indeed, for N <C L one finds from (|3.11|) that £ ~ 1/L and the condition 
£ < 1 is satisfied. For N ^> L one gets from (j3.11|) 

1 N 

£~-ln-. (3.15) 

One recognizes that similar parameter £ str naturally arises on the string side and controls there 
the transition between two different regimes in (|1.4jl . For £ ^> 1 one has to take into account 
that the two cuts collide at the origin and, as a consequence, the analytical properties of S'(x) 
are different. The simplest way to see this is to examine the Taylor expansion of the leading term 
Sq(x), Eq. (|3.6|) . around the origin and retain the terms most singular for b — > 



S o( x ) ~ ^ 



1 



x 2 /b 2 



JL 

2b 2 



3 x 2 5 x 4 
1 + 4 K 2 + 8 ¥ + 



(3.16) 



As follows from ()2.14j) and ()2.33)1 . higher order terms of this expansion contribute to subleading 
corrections to the energy e = i f + dx S"(x). It is easy to see that the series for i J x + dx Sq(x) 
runs in powers of £ = and it is only convergent for £ < 1. 

Based on our analysis we expect that the semiclassical expansion of the energy ()2.33)1 will be 
divergent for £ ^> 1. Let us calculate the first subleading correction to the energy ()2.33|) defined 
by S"(0). The function S[(x) satisfies the relation ()2.18|) which involves yet unknown function 
Ti(x). To determine S[(x) one requires that it should have the same analytical properties as 
the leading term of the semiclassical expansion S' Q (x), i.e., to be an analytical function in the 
complex plane with the cuts [—a, —b] U [b, a]. Discontinuity of S[(x) across the cut defines the 
correction to the distribution density of the Bethe roots [2H]. Since tq(x) and Ti(x) are entire 
functions while p'(x) and smp(x) change sign across the cut, sinp(x + zO) = — smp(x — iO), the 
function S[(x) satisfies the relation 

S[(x + zO) + S[(x- iO) = -p'(x + z'O) cotp(x + z'O) , (3.17) 

with x G [—a, —b] U [b, a]. In addition, its asymptotic behavior at the origin and infinity is fixed 
by Eqs. (gilBD and 

^(0) = 0, S[(x) x ~°°±. (3.18) 

x 6 

The solution to the resulting Riemann-Hilbert problem looks like (c.f. j37j ) 



x 



v/(x 2 -a 2 )(x 2 -b 2 ) J b vr 



dz y/{a 2 - z 2 )(z 2 -b 2 ) 



= n 5 „ n ^? / ~ V o_ \ L v'{z + zO) cotp(z + zO) , (3.19) 



x* — 
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where £ +iO' indicates that the momentum is evaluated on the upper sheet of the Riemann surface. 
Replacing p'(z + iO) by its expression ()3.2j) . we obtain the following representation for the first 
subleading correction to (|2.33|) 



Ei 



2 JL 

L 



S?(0) 



2/? 
nabL 



r dz 




L * 






z 2 



coth(ip(z + zO)) . 



(3.20) 



For our purposes we would like to determine the asymptotic behavior of E\ for N ^> L, or 
equivalently a — » l/(2m), b — » and (3 — > 0, Eq. (|3.13|) . Changing the integration variable to 
z — > oz, one finds in the limit 6^0 



El 



2/9 


r°° dz 


" /3a 1 " 


Tcab 2 L 


A ^ 


6 z 2 



coth.{ip{bz + iO)) 



Ail 



+ . . . . 



(3.21) 



Here in the second relation we took into account that for (3 — ► the cuts collide at the origin 
and the momentum scales as p(z + iO) ~ 2imlnz for z — > 0, Eq. 1)2.47)1 . Comparing ()3.21jl with 
a similar relation for the leading term Eq ~ ft 2 /{Lb 2 ), Eq. (|3.7)1 . one finds that 



= 0(17/6) = O(0 



(3.22) 



where ry/6 = 1/(6JV) = (2m/s) \n(N/L)/L in the limit > L, Eqs. ()2~T0)) and 

We conclude that foiN^L the semiclassical expansion of the energy, e = Eo + ^£i + 0(r^ 2 ), 
Eq. ()2.33)) runs in powers of £ and it is only convergent for £ < 1. Together with ()2.2j) this 
leads to the expression for the one-loop anomalous dimension, Eqs. f)1.7p and (|1.8|) . with c x (1 = 
—2m/ (3ns). 



3.3 Beyond semiclassical expansion 

We demonstrated in the previous section that the semiclassical expansion breaks down for £ ^> 1 . 
The reason for this is that the cuts collide in this limit and the semiclassical expression for S(x) 
ceases to be analytical at x — 0. Substituting /3 = into ()3.6|) . one finds (for m — 1) 

S£(s) = = (3.23) 

and, therefore, S' (x) ~ =F2ilnx for i -> on the upper and lower sheets, respectively. It is easy 
to see that the corresponding transfer matrix t (x) = 2cosp(x) is given by (|2.44j) . We remind 
that tq(x) plays the role of the potential in the Baxter equation (|2.3jl . The fact that tq(x) is 
singular at x = implies that the semiclassical expansion (|2.15j) breaks down at the origin. In this 
section, we present an approach which allows one to construct asymptotic solution to the Baxter 
equation ()2.3|) and evaluate the energy for N > L. It is complementary to the semiclassical 
approach and takes a full advantage of the above mentioned singularity of the transfer matrix 
Tq(x) at x = 0. A detailed account on this approach can be found in Ref. [To] . 
To begin with, we rewrite the Baxter equation ()2.3j) as 

(u + is) L (j>{u) + ^ ~ ' g) . L = t L (u) (3.24) 
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where the notation was introduced for 

Q(u + i) 

<PW = nl v ■ 3.25 

To evaluate the energy ()2.8|) one has to construct solutions to ()3.24|) for u ~ ±is. Notice that 
the dressing factors {u ± is) L in ()3.24j) vanish for u — > =pis indicating that one of the terms in 
the left-hand side of (13.24)1 becomes anomalously small and can be neglected. It turns out that 
in the thermodynamic limit, L ^> 1, the same approximation can be performed not only in the 
vicinity of u = Us but in the whole region u = O ((N + Ls)°). To see this we note that ti{u) is 
given by ()2.4|) with the charges that scale as ~ (N + Ls) k . This suggests that ^ 1 f° r 

u = O ((N + Ls)°). Indeed, in the semiclassical approach, for u = x/r), one finds from (j2.14|) . 
(j2.17|) and ()3.23|) that the transfer matrix takes the form 

t L (x/rj) = 2(x/r]) L cosp(x) + 0(rj) . (3.26) 

For u = O(L ), or equivalently x ~ rj, the momentum can be replaced by its asymptotic behavior 
at the origin: p(x) ~ (3/x for £ <^ 1 and p(x) ~ 2i\nx for £ ^> 1 leading to |t^(x/?7)| 3> 1. 

Let us return to Eq. (|3.24j) and take into account that tj^iu) takes large values for u = 
O ((N + Ls)°) both for £ < 1 and £ 3> 1. Requiring the left-hand side of ()3.24|) to be as large 
as ti(u) one finds that either <p{u) ^> 1, or <p(u -i) <C 1. In both cases, one of the terms in the 
left-hand side of ()3.24|) can be safely neglected and one arrives at two different equations 

(u + is) L Q + (u + i) = t L {u)Q + {u) , (3.27) 
(u - is) L Q_(u -%) = t L {u)Q^{u) . 

Having solved this system, one can construct an asymptotic solution to the Baxter equation (|2.3|) 
as a linear combination of Q+{u) and Q-(u) 

QW(u) = c + Q + {u) + c_Q_(m) . (3.28) 

We would like to stress that Q^' (u) does not satisfy the Baxter equation (|2.6|) . but asymptotically 
approaches its solution Q{u) in the region u = O ((N + Ls)°). Equations ()3.27|) were obtained 
under the assumption that {u + is) L Q + {u + i) ^> (u — is) L Q + {u — i) and {u + is) L Q-{u + i) 
(u — is) L Q_(u — i), respectively. Together with (|3.27|) it can be expressed as the following relation 
for the transfer matrix 

t L (u + |) t L (u - |) > \u 2 + (i - s) 2 ] L . (3.29) 

For u = O ((N + Ls)°) it is equivalent to \t L (u) \ > 1. 

To solve (|3.27|) one introduces into consideration the roots of the transfer matrix ()2.4|) 

L 

t L {u) = 2\[{u-8 n ). (3.30) 

71=1 

It is known that for polynomial solutions to the Baxter equation, Eq. ([2.6)1 . the roots take real 
values only, Im 5k = |2E] • Matching ()3.30)) into ()2.4j) one finds that they satisfy the sum rules 

X> = o, X>' = -^ 2 ' II 6 - = ■ ( 3 - 31 ) 

n=l n=l n=l 
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For even solutions to the Baxter equation, Q{u) = Q(—u), the roots appear in pairs 5 n = 
—Sl-u+i- Making use of ()3.30|) . it is straightforward to verify that the solutions to ()3.27|) are 
given by 



, T(—iu + s) " 1J - T(iu + s] 

n=l n=l v ' 



To fix the constants c± in ()3.28|) one examines the relation for the quasimomentum in Eq. (|2.8j) 
and substitutes Q{u) by its asymptotic expression (|3.28p . Taking into account that (for g — > 0) 

Q+(-is + ~ q l , Q^{is-g)~Q L (3.33) 

one finds from (£BJ) and (13~251) 

Q (as) (^) = c+ g+gg) ( „ ^ 

Q(»)(-is) c_Q_(-is) ' 1 ' 



e i8 



Therefore, for cyclically symmetric states e te = 1 the asymptotic solution to the Baxter equation 
is given up to an overall normalization factor by 

QM( U ) = Q + ( u )Q_(-i s ) + Q_( u )Q + (is) . (3.35) 

We remind that this relation is only valid in the region u = O {{N + Ls)°). In a similar manner, 
one uses (I3.33J) to evaluate the energy ()2.8|) as 

e {as) = i (in {is))' -i (In {-is))' = i(lnQ+(is))' - i(lnQ_(-is))' . (3.36) 

Replacing Q±(u) by their expressions ()3.32|) one finds the following remarkable expression for 
the energy in terms of the roots of the transfer matrix ()3.3())1 [2EI E3 E] 

L 

e (as) = 21n2 + ^[^(s + ^ n )+^(s-^„)-2^(2s)] , (3.37) 

n=l 

where ^(x) = <ilnr(x)/c?x. The interpretation of (J3.37j) in terms of classical SL{2) spin magnet 
and properties of the Wilson operators (|2.1|) in gauge theory can be found in Ref. 

Since the roots 5 n are functions of the conserved charges, Eq. (j3.31|) . the relation (j3.37|) 
establishes the dependence of the energy on q 2 , . . . , qi- To check ([3.37)1 . we compared the asymp- 
totic "dispersion curve" e^{q2, ■ ■ ■ , q£) with the exact one e^{q 2 , ■ ■ ■ , <?l) coming from the 
Bethe Ansatz solution (|2.9|) . We found that for L — 10 and s = 1/2 the accuracy of (|3.37|) . 
e (as)/ e (ex) _ 1; i ncr eases from -4.6 x 10~ 5 for A" = 2, to -2.6 x 10~ 6 for A^ = 10 and to 1.9 x 10~ 10 
for A" = 100. We conclude that Eq. ()3.37|) describes the exact eigenspectrum with a high accu- 
racy throughout the whole interval of A" including the region A" ~ L in which the semiclassical 
approach is applicable. Indeed, it is straightforward to show that the relation ()3.37|) coincides 
with the semiclassical expression (J2.33j) for N ~ L (see Appendix). 

The charges 93, . . . , ?l have to satisfy quantization conditions. In the method of Baxter 
Q— operator they follow from the requirement for Q{u) to be polynomial solutions to the Baxter 
equation, Eqs. (|2.3|) and ()2.6|) . The roots of Q{u) scale in the thermodynamic limit as u = \j ~ L 
and, therefore, they lie outside the applicability range of ()3.35|) . This does not allow us to impose 



21 



the polynomiality condition for Q^(u). We shall require instead that Q^(u), given by ()3.35j) . 
should be regular on the real u— axis inside the region u = O ((N + Ls)°). This property is not 
warranted since both Q+(u) and Q-(u) develop poles on the real u— axis originating from the 
product of T— functions in the numerator in the right-hand side of Eq. ()3.32|) . The poles are 
located at zeros of the transfer matrix f)3.30|) . u = 8j with j — 1, . . . , L. Requiring for Q( as )(u) to 
have zero residues at u = 8 n with 8 n = O ((N + Ls)°), one gets from ()3.35|) and ()3.32|) 



-2is n tt T(-i5 n + i5 



T(s — iS r , 



L L 



T(s + i8i 



V(i8 n -i8 j ) [T(s + t5 n )\ l = lr(a-%)' 



The total number of roots of the transfer matrix ()3.30|) equals the length of the spin chain L. 
As we will see in a moment, in the thermodynamic limit L ^> 1 and iV ^> L all roots can be 
separated into two different groups depending on their scaling: "small" roots 8 n = O(N ) and 
"large" roots 8 n = O(N). It is important to realize that the quantization conditions ()3.38|) should 
only hold for the small roots while the product over j entering both sides of (j3.38|) involves all 
roots. 

The total number of small roots depends on the value of the parameter (3 = Ls/(N + Ls). 
For < (3 < 1, that is, within the applicability range of the semiclassical expansion, the integrals 
of motion scale as in Eq. (|3.1U|) . Together with (|3.31|) this implies that the roots behave as 
5 n = O(N) and, therefore, there are no small roots. For f3 — > higher charges g^, qL-i, ■ ■ ■ take 
anomalously small values (j3.1U|) indicating that the small roots are there. Let us demonstrate 
that their total number equals L — 2m with positive integer m entering the right-hand side of 
(J3.3j) . One makes use of ()3.2fi|) to deduce that the large roots of the transfer matrix satisfy 
cosp(5jr]) = + 0(f)). Recall that the momentum verifies (cosp(x)) 2 > 1 on the cuts (for 
b 2 < x 2 < a 2 ) and takes the values cosp(x) = ±1 at the double points x = X2j, Eq. ()3.9|) . Since 
cosp(x2j) = — cosp(x2j-2), the roots of the transfer matrix ()3.26|) should lie on the real axis in 
between the branching points, X2k < djV < %2j-2, away from the cuts [—a, —b] U [b, a] (see Fig.EJ). 
Going over to the limit (3 = Ls/(N + Ls) — > 0, or equivalently N ^> L, one takes into account 
that the edges of the cuts scale as a = l/(2m) and b = f3/(— 2m In/?) — > 0, Eq. (|3.13|) . As a 
result, all roots of the transfer matrix 5ji] ~ Sj/N "trapped" inside the interval [—b, b] are to be 
small while the roots belonging to the intervals (— oo, —a] U [a, oo) are to be large 

| 5 (large)| > l^small), < J_ ^ ^ 

2m 4m£ 

with £ = \n(N/L)/L. Thus, the total number of small roots, L — 2m, equals the number of 
double points on the interval [—6, b] and, as a consequence, the total number of the large roots 
is 2m. 

Let us consider the minimal energy state with m = 1. For N ^> L, the transfer matrix has 
two large roots 5i = —5l and L — 2 small roots 8k = —8L-k+i with k — 2, . . . , L — 1. It follows 
from (j3.31|) that the large root is given by 

8, = (- q2 /2) 1/2 [1 + 0(1/L)) ~ ^±=^ (3.40) 

with q 2 defined in f!2.5|) . The small roots satisfy the quantization conditions (j3.38J) . Separating 
in ()3.38|) the contribution of the large roots with j = 1 and j = L, one can rewrite ()3.38|) in 
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terms of the small roots only (for n — 2, . . . , L — 1) 



T(s - i5 n ) 
T(s + i5 n ) 



L L-l 



„T(l + i6 n -i6 s ) 

11 r(i + 1 } 



Here we took into account that the product over j in the right-hand side of 1)3.38)1 equals 1 due 
to pairing of the roots 5 n = —di-n+i- Taking the logarithm of both sides in 1)3.41)1 . one rewrites 
the quantization conditions as 

L-l 

S n \n(-q 2 ) + LargT(s - i5 n ) + ^argT(l + i5 n - iSj) = -k n (3.42) 

i=a 

where integers k 2 > k% > . . . count the branches of the logarithms and satisfy k n = — k,L- n +\. 
In addition, they have a parity opposite to that of L, that is k n = even/odd for L = odd/even. 
Notice that in distinction with the Bethe ansatz equations ()2.7j) . the number of relations in ()3.42j) 
does not depend on the total spin iV but only on the length of the spin chain L. 

To evaluate the energy ()3.37)) one has to solve the system ()3.42j) and determine the set of 
small roots. The resulting expression for e^ s \ Eq. f[3~37|l . depends on integers k n and as we show 
below (see Eq. (J3.48)) ) it takes minimal value for the occupation numbers k n = L + 1 — 2n. A 
systematic analysis of the quantization conditions ()3.42jl will be given elsewhere. The system of 
equations ()3.42j) can be easily solved in two limits £ ^> 1 and £ < 1. 

For £ ^> 1 (or ln(iV/L) ^> L) one deduces from ()3.39|) that solutions to ()3.42|) have to satisfy 
\8 n \ -C 1. In this case, one expands the T— functions in ()3.42|) in powers of <5's and one finds after 
some algebra 

7ck n /2 

6n = H-q 2 ) + (L - 2)^(1) - ms) + • • • (3 - 43) 

with n — 2, . . . , L — 1. One verifies a posteriori that this relation is in agreement with ()3.39|) . 

For £ < 1 (or 1 < \n(N/L) < L), one obtains from (l3~39l that |<J n | < l/(2f) and, therefore, 
not all roots verify the relation \S n \ 1. Still, the relation ()3.43|) is valid for the roots with small 
absolute value \5 n \ -C 1 (with n ~ L/2). For roots \8 n \ ~ l/(2£) (with n ~ L°) one expands the 
T— functions in (|3.42j) in inverse powers of 5's and obtains 

5 n \n(-q 2 )-5 n \n5 2 n + ... = ^k n , (3.44) 

with k n = O(L). From this relation one obtains the relation 5 n ~ 7rA; n /(41n(A^//c n )) which is in 
agreement with (|3.39|) . 

To test the quantization conditions ()3.42j) we compared solutions to ()3.42j) for k n = L + 1 — 2n 
with the exact, Bethe ansatz expressions for the roots of the transfer matrix corresponding to 
the minimal energy state with the quantum numbers s — |, L — 10 and iV = 100 and observed 
a good agreement (see Table HJ). In agreement with (j3.43J) and ()3.44|) . the small roots 5 n vary 
linearly with n close to the origin and deviate from the linear behavior close to the end points. 
Moreover, for £ 3> 1 all small roots scale linearly with n in agreement with (J3.43|) (see Fig. 

In a similar manner, we solved quantization conditions (J3.42)) for s — |, L — 10 and < iV < 
100 and compared the resulting expression for the energy (|3.37j) with the exact expression (|2.9j) as 
shown in Table El We found that for iV > L the asymptotic expression ()3.37|) approximates the 
exact result with an accuracy better than 0.002% while the semiclassical approach significantly 
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Si 


5 2 


Ss 


s± 


s 5 


E 


exact 


73.897289 


0.5297596 


0.3443964 


0.1937491 


0.06253897 


7.3790455 


asym 


73.900271 


0.5297573 


0.3443955 


0.1937487 


0.06253886 


7.3791719 



Table 1: Comparison of the exact roots, 5 n = — Sl-u+i, an d the energy, E, for s = 1/2, L = 10 
and N = 100 with the asymptotic expressions obtained from Eqs. (|3.41|) and (J3.37)) . respectively. 



N 





10 


20 


30 


40 


50 


60 


70 


80 


90 


exact 





2.2766 


3.5069 


4.3644 


5.0266 


5.5678 


6.0262 


6.4243 


6.7765 


7.0924 


asym 


0.0907 


2.2851 


3.5097 


4.3657 


5.0274 


5.5683 


6.0265 


6.4246 


6.7767 


7.0926 



Table 2: Comparison of the exact energy, e, with the asymptotic expression obtained from 
Eqs. fHH) and (13371) for s = §, L = 10 and different total spin N. 

overestimates the value of the energy (see Fig. 0J). Remember that the minimal energy states 
carry even Lorentz spin only, Eq. (|2.31|) . 

Let us examine the expression for the energy ()3.37j) in the limit £ ^> 1, or \n(N/L) ^> L. We 
remind that the semiclassical expansion breaks down in this region. Equation ()3.37|) involves the 
sum over the roots of the transfer matrix ()3.30|) . As before, we separate them into two groups. 
According to f)3.39|) . the large and small roots scale as Sj ~ N and Sj ~ l/£, respectively. This 
allows one to replace ip(s ± iSj) in ()3.37|) by its asymptotic behavior at infinity and at the origin, 
respectively. In this way we obtain 

e = 21n2-2L^(2s)+2^1n|5 j | + [ 2 *P( S ) ~ *P"( S ) tf] +■•■ , (3.45) 

large small 

where the ellipsis denotes subleading terms. The number of large roots equals 2m, Sj = —Sz-j+i 
(with j — 1, . . . , m) and, therefore, the leading asymptotic behavior of the energy for £ ^> 1 is jS] 

e = A\n(5 1 5 2 ...5 m ) + O(N ) = 2\n\q 2m \ + O(N ) , (3.46) 

where we took into account the relations (|3.31|) and q2 m = q2mN 2m was defined in (I2.48J) . In this 
way we obtain the relation 

£ = 4mlniV + e>(iV ), (3.47) 

which is in an agreement with (J2.42|) . The minimal energy corresponds to m = 1. 

For m = l one has just two large roots 5i = — 5l and L — 2 small roots 5 n = —5L~ n +i (with 
n = 2, . . . , L — 1) defined in Eqs. ()3.40|) and ()3.43|) . respectively. Therefore, for £ ^> 1 one gets 

e = 4 In N + 2L [^(s) - ^(2s)] + !>' + ••• (3-48) 

71 = 2 

This expression depends on the set of integers k 2 > k 3 > . . . (with k n = —kL-n+i) defined in 
1)3.42)1 . Since ip"{s) < for s > 0, the minimal value of e corresponds to k n = L + 1 — 2n 

e = 4 In N + 2L Ufj(s) - M2s)] + 1 - 9 l '' +... (3.49) 
L v 481n 2 A^ v 1 
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Figure 4: Left panel: "small" roots of the transfer matrix S n (with n — 2, . . . , L — 1) for s — ^, 
L — 10 and two values of the spin iV = 10 2 and N = 10 10 . Crosses stand for the solutions to 
the quantization condition ()3.42|) and the lines correspond to ()3.43|) for k n = L + 1 — 2n. The 
roots tend to approach the line as £ = \n(N/L)/L increases from £ = 0.23 to £ = 2.07. Right 
panel: the minimal energy for s = |, L = 10 and the total spin < N < 100. Crosses denote the 
exact values, Eq. (j2.9|) . while the solid line stands for the semiclassical expression f)3 . 7j) . We do 
not display the data for the asymptotic energy ()3.37|) with roots deduced from the quantization 
conditions ()3.42|) . since they are not distinguishable from the exact spectrum. 



This relation defines the minimal energy for £ ^> 1. The energy of the excited states is described 
by (|3.48|) with another set of k— integers but of the same parity. The lowest lying excited state 
has the same fc's as the minimal energy state except &2 = —ki-i = L — 1. It is separated from 
the latter by the distance 

2 In N ln 2 iV V ' 

This relation defines the level spacing in the spectrum of anomalous dimension of the Wilson 
operators (|2.1j) close to the minimal anomalous dimension trajectory. 

So far our discussion was limited to one loop. Remarkably enough the logarithmic behavior 
of the anomalous dimension ()2.2|) persists to all orders of perturbation theory. Higher order 
corrections to ()2.2)1 merely modify the coefficient in front of In N replacing A by an infinite series 
in the coupling constant, 

7 (A) = A [ 4 In N + 0(N )} + 0(\ 2 ) = 2T cusp {X) In JV + O(N ) , (3.51) 

with r cusp (A) being the cusp anomalous dimension. This relation holds true in a generic Yang- 
Mills theory for arbitrary values of the coupling constant A including the strong coupling regime. 
In the latter case, one can apply the gauge/string correspondence to obtain a prediction for the 
cusp anomalous dimension at strong coupling in the M = 4 SYM theory. In the next section, we 
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shall trace the origin of logarithmic scaling of anomalous dimension in the dual picture of folded 
string rotating on the AdS space. 

4 Classical strings in ACIS3XS 1 

Let us turn to the analysis of anomalous dimensions of composite operators (j2.1j) on the string 
side. According to the gauge/string correspondence the Wilson operators are mapped into certain 
string states whose energy is identified with the scaling dimensions of the former. It is known that 
in the M = 4 SYM theory, Wilson operators with the minimal anomalous dimension discussed 
in the preceding sections are dual to a single-folded string rotating in the AdSsxS 1 sector of the 
target space of the type IIB string theory [HUH]- 

4.1 Folded rotating string 

For the folded closed string spinning both in AdS3 and S 1 , the relevant bosonic part of the 
superstring action reads 

S = ~ / d 2 i V- det \\G MN {X) d a X^d h X^\\ , (4.1) 

where the 't Hooft coupling A = g 2 N c is related to the product of the radius of anti-de Sitter 
space R and the string tension 1/a', y/X = R 2 /a', and the metric in the target space looks like 

ds 2 = G M N<iX M dX N = - cosh 2 p dt 2 + sinh 2 p d6 2 + dp 2 + dip 2 . (4.2) 

Here t, p are 9 are the global time, the radial coordinate and the angle, respectively, on the 
AdS3 space and <p is the angle on a large circle of S 5 . The worldsheet coordinates of the string 
£ a = Z, 1 ) are chosen in such a way that 

t = ?, 9 = ue + Z\ <p = v?. (4.3) 

Here u is the angular velocity on AdS and v is a boost parameter of the center-of-mass of the 
folded string on S 5 . For the rigid folded string the radial variable p does not depend on £° 
and is periodic in that is, p(£ x + 2ti) = p^ 1 ). The dependence p = p(£ x ) is determined 
by the classical string equations of motion [UJ HTJ EE!- Their solution describes a folded string 
configuration which is sewed out of four segments running along the radial direction from the 
center of AdS space, p = 0, towards its boundary by a distance po, which depends on the angular 
velocities entering (|4.3j) 

9 to 2 — v 2 

coth 2 p = ^ r >l. (4.4) 

The string has two spikes located at p = po which are responsible for the logarithmic scaling of 
the anomalous dimension of Wilson operators at strong coupling |%| I38j. 

Within the gauge/string correspondence, the anomalous dimension of the operators (J2.1)) is 
related to the energy E of the classical rotating string as 

E = iV + L + 7 (A), (4.5) 
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where the Lorentz spin N and twist L are translated into the angular momenta on the AdS3 and 
S 1 spaces, respectively. From the Nambu-Goto action (j4.1|) . one finds these charges as 

7T Vl-r 7T y/l — IS 2 ^ y/l — V 2 

(4-6) 

where K(x) and E(x) are the elliptic functions of the first and second kind, respectively, and the 
auxiliary (negative valued) parameter x is related to the distance po from the center of anti-de 
Sitter space to the end of the string extending to its boundary, 

X = -sinh 2 p - (4.7) 
Excluding uj and v in favor of Xi Eqs. f!4.6|) can be cast into the parametric form 



Using these relations, one can analyze the anomalous dimension r y = E — N — LaX strong coupling 
in the semiclassical limit of large angular momenta L, N — > oo with N/ y/X and L/yfX kept fixed. 
In this way one finds from (|4.8|) that the two limiting cases N <C L and N ^> L correspond to 
the short and long strings, po 1 and p ^> 1, respectively. 

It is well-known jHHH] that in the short string limit, po "C 1, or equivalently (— x) <S 1, the 
anomalous dimension exhibits the BMN scaling 

7 (A) = L [a' 7 (0) + (A') 2 7 (1) + • • .1 , (4.9) 

with 7^°' ) , 7W, . . . being functions of the ratio N/L and A' = X/(irL) 2 determining the BMN 
coupling constant. Substituting x — Xo + ^'Xi + • • • into ()4.8|) and matching the coefficients in 
the expansion of both sides of (|4.8jl in powers of A', one finds for N/L <C 1 [TU1 ITT| ITHj 

7«°» = YT + ..., 7 « = - TT + .... (4.10) 

The lowest order term in ()4.9|) matches the one-loop result on the gauge theory side, Eq. ()3.12|) . 
for s — | and m = 1 

7 (0) = ^- (4-11) 

In a similar manner, in the long string limit, po — * 00 ; or equivalently (— x) — > 00, assuming 
that x has a regular expansion in powers of A', one finds from Eq. (|4.8jl that the leading order 
parameter xo scales for N/L ^> 1 as 

N N , . 

Xo = -^ln r . (4.12) 

This leads to the following expressions for the coefficients of the BMN expansion of the anomalous 
dimension (JO)) for N/L > 1 

7 (0) = |ln 2 (- Xo ) + . . . = ±ln 2 (iV/L) + . . . , (4.13) 
7 (D = _ 1 i n 4 ( _ Xo) + ... = _! ln 4 (iV/L) + . . . , 
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where ellipses denote terms subleading for (— xo) oo. In agreement with (j4.11j) . the relation 
(14.13)1 matches the one-loop expression for the anomalous dimension on the gauge theory side, 
Eq. (|3.14|) . Moreover, one can show 18, that the relation (|4.11j) (upon the Landen transformation 
defined below in Eq. (14.33)) ) holds for an arbitrary value of N/L with 7^ given by ()4.8|) and 
(@~HD and e defined in fUJ) and (HOI) . 

So far we see no trace of expected scaling of the anomalous dimension 7(A) ~ In N for N ^> L. 
We recall that on the gauge theory side the semiclassical expansion for the one-loop energy e was 
divergent for N ^> L and in order to recover the logarithmic scaling e ~ IniV we had to resum 
the entire semiclassical series in powers of the parameter £ = hx{N / L) / L. It turns out that a 
similar phenomenon happens for the series ()4.9|) on the string side although the parameter of the 
semiclassical expansion is different and equals 



str 



A 

1? 



In 2 



N 



(4.14) 



It can be easily identified by comparing the contribution to (J4.9|) from 7^ and 7W, Eq. (|4.13|) . 

To sum up the infinite series in ()4.9|) for £ str > 1 we return to ()4.8|) and rewrite the first 
relation as 

VT- 4A' X K 2 (x)-1 + A 7 



7(A) = L 



with A' = A/(ttL) 2 and 



A7 



4A'K(x)[E(x)-K( X )] 



y/1 - AX'x^(x) + y/l- 4A'(* - 1W( X ) 



(4.15) 



(4.16) 



In the limit of a long string, for (— x) —* 00, the expression for 7(A) can be expanded in powers of 
4A'(— x)K 2 (x) ~ A' In (— x)- Examining the expression for A7 one finds that its contribution to 
7(A) is suppressed by the factor \E(x) ~^(x)}/ (~ X^-(x)) ~ 1/ ln(— x) compared to the first term 
in the square brackets in ()4.15|) . This implies that in the limit (— x) 00 with A' ln 2 (— x) = fixed 
the leading asymptotic behavior of ()4.15|) reads 



7(A) 



1 + A'ln 2 (- X )-1 



(4.17) 



This relation resums all double-logarithmic corrections ~ L[X' ln 2 (— x)] n to the anomalous di- 
mension 7(A) to all orders n at strong coupling. In particular, for n = 1 and n = 2 it reproduces 
(I4.13)) . For A'ln 2 (— x) <^ 1 one expands the square-root in ()4.17|) and arrives at a BMN-like 
expansion (j4.9j) . At the same time, for A'ln 2 (— x) ^> 1 the relation ()4.17|1 leads to the expression 



7 (A)=WA'ln(-x) + 



VA 



7T 



H-x) + 



(4.18) 



For (— x) —> 00, the dependence of x on the ratio N/L and the BMN coupling A' follows from 
the second relation in ()4.8|) 



4x 2 



17_ 

N 2 



+ A' 



(4.19) 



For £ srt < 1 and N > L one finds from ()4.19j) and (j4.14|) that x 



Xo 



"2Z m T an< ^' therefore, 



7(A) = L 



1 + 



A 



In 2 {N/L) - 1 



(4.20) 
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For £ srt > 1 one finds from P~T9jl and (jHHl) that (-%) ~ N/{2L\f\') ~ iV/v 7 ! and, therefore, 



7 (A) = — ln(iV/VA) + ... (4.21) 

7T 

This relation reproduces the correct asymptotic behavior of the anomalous dimension in the 
regime ()1.4|) and it is in agreement with the results of Refs. |d1 ITU I3S]. 

The coefficient in front of IniV in ()4.21|) determines the leading asymptotics of the cusp 
anomalous dimension in the strong coupling regime in the M = 4 SYM theory, Eqs. and 
(jl.2j) . It follows from the relations (J4.2Uj) and (|4.21j) that quite remarkably the cusp anomaly at 
strong coupling can be obtained by resumming double logarithmic terms ~ L[Aln 2 (A^/L)/L 2 ] fc 
in the expansion of anomalous dimensions of operators of higher twist L and large Lorentz spin 



4.2 Two-cut solution in string sigma model 

In the previous section, we observed that for N ^> L the dependence of the anomalous dimension 
7(A) on the coupling constant is different for £ str < 1 and £ str 3> 1. In the former case, 7(A) 
has a BMN-like expansion in powers of A/L 2 , while in the latter case 7(A) is not analytical in 
A and scales logarithmically ~ v^A In N. The reason for this non-analyticity is that for £ str ;§> 1 
the end-points of the rotating string approach the boundary of the AdS space and the dominant 
contribution to the energy of the string, or equivalently the anomalous dimension ()4.5j) . comes 
from the vicinity of two spikes, 7 (A) ~ (y/X/ir) ■ 2p . Here the radial coordinate of the spikes 
scales for £ str > 1 as p ~ ^ ln( — x) ~ | ln(iV/A/A), Eqs. ()4.7|) and (j4.19|) . The phenomenon is 
rather general [HI EE] and it holds true for classical string configurations with an arbitrary number 
of spikes n. In that case, each spike provides a logarithmic contribution to the energy and the 
anomalous dimension scales for £ str > 1 as 7(A) = ny/\/(27r) \n(N/\^X). At first glance, this 
mechanism is rather different from the one in gauge theory. We recall that in gauge theory, to 
one-loop order, the logarithmic scaling of anomalous dimension for N ^> L arises due to collision 
of cuts for the spectral curve of the spin chain, Eq. (j2.19|) . It is known that the classical equations 
of motion for the string on the AdSsxS 5 background are completely integrable [221 HI] and their 
solutions are parameterized by the spectral curves. Moreover, for the strings on the AdS3xS 1 
part of the target space the spectral curve can be identified as a complex hyperelliptic curve 
For the folded rotating closed string discussed in the previous section, it is given by the elliptic 
curve with symmetric branching points on the real axis [21] 

r str : y 2 = (x 2 -a 2 stI )(x 2 -b 2 stI ), (4.22) 

with a s tr and 6 str taking positive values, 6 s t r < a s tr- In this section, we shall translate different 
asymptotic behavior of the anomalous dimensions for N ^> L, Eqs. (|4.2U|) and (J4.21j) . into 
properties of the curve (|4.22j) and reveal the mechanism responsible for the logarithmic scaling 

Similar to the classical SL(2) spin chain, the classical string equations of motion admit the Lax 
representation and they can be solved exactly by constructing the Baker- Akhiezer function [231 
[22]- As before, the Bloch-Floquet multiplier for this function gives rise to the (quasi) momentum 
p(x) which is the generating function for the conserved charge including the energy. For the 
folded rotating string configuration, p'(x) is an analytical function in the complex plane with 
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two symmetric cuts [— a str , — b stT ] U [b s tr,a stT }. It is uniquely fixed by the requirement that dp = 
p'(x)dx should be a meromorphic differential on the complex curve (|4.22jl satisfying the following 
conditions on the upper sheet of r str [21] : 



Prescribed asymptotics at infinity and at the origin 

dxE + N 



dp 



x— >oo 
~ — Z- 



x- 



Single-valuedness condition 



dp = 



f'str 



Double poles at x — ±y/X/(irL) = vA 

dp ~ dx 



, 7 2 E — N 

dp ~ —dx— — 

A .L 



-7rm . 



x±VA 



^2 



+ 0((x ± v/A 7 ) ) 



(4.23) 



(4.24) 



(4.25) 



The resulting expression for the differential dp takes the form 

y+ v+ y+ 



, dx 
dp = — 

y 



(x — vA 



M2 



X 



A' (x + VA') 2 x + V A 



C 



(4.26) 



where y = y(x) is defined in ()4.22j) . y + = y(vA') and = y ; (VA 7 ). 

Equation (|4.26j) depends on three parameters a stT , 6 str and C. They are fixed by the normal- 
ization conditions (J4.23|) and (|4.24|) as 



1 



-'str 



c 



mK(r) 
ma str 



fl 2 " 

u str 



E(r) - 4~ K ( r ) 



u str 



-1/2 



(4.27) 



with the modular parameter r = 1 — 6 2 tr /a 2 tr . In addition, one finds from ()4.23j) the following 
expressions for the ratio N/L and for the anomalous dimension 7(A) = E — N — L 



N/L 



m 


E(r) ^a str + 










-V 


- K(r) (&str H 




~ ~2 




bstrj 




a str / . 






n A'" 






= m 


K(r) 6 str - E 




-1, 








°str_ 







(4.28) 



with 6 str defined in ()4.27|) and a str = festr/ Vl ~~ r - 

Let us examine the dependence of the anomalous dimension on the coupling constant A' = 
A/ (ttL) 2 . Assuming that fe str and r both admit a regular expansion in powers of A', one substitutes 
into (1071) and KM 



''str 



&str + A' £"str 



+ . . . 



r = r 



(°) + A / r« + ... 



(4.29) 
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and matches the coefficients in front of powers of A'. In this way, one obtains to leading order 



6 (o) _ 1 N _ 1 



y str 



ml(r(°)) ' L 2 - tWK(t(°)) 



(4.30) 



and all subleading corrections to (J4.29)l are expressed in terms of the leading terms. Then, the 
first few coefficients of the BMN series for anomalous dimension ()4.9|) are given by 



9 

m 



7 (o) = ^L K ( r (o)) [(2 - t< >)K(t(°>) - 2E(r(°))] , (4.31) 
7 (D = ^K 3 (r(°)) [(4(2 - r<Vl - r(°) - (r< f )l(r(°») - 8^1 - r(»)E(r(°))] . (4.32) 

Together with the second relation in 1)4.30)1 . they define the parametric dependence of the anoma- 
lous dimension ()4.9j) on N/L. 

Equations ()4.30|) determine perturbative corrections to the anomalous dimension of long scalar 
operators in the M = 4 SYM theory. According to 1)4.11)1 . it should match at one- loop order a 
similar asymptotic expression ()3.7jl obtained on the gauge theory side within the semiclassical 
approach. The conformal spin of scalars equals s — 1/2 and the scaling parameter j3 = L/(L + 
2N), Eq. (j2~T0|) . is given by (3 = Vl - r(°) K(t (0) )/E(t^) in agreement with (Q. Then, one 
observes that 7^°- ) and £ given by (j4.31)l and ()3.7)) . respectively, verify the relation (j4.11j) . The 
expressions for the branching points b^l and b, defined in ()4.30)) and ()3.4jl . respectively, are 
different, = /3/2, but the agreement can be restored through the rescaling of the local 

complex parameter x in the definition of the curve (j4.22j) . x — > x/3/2. 

The functional form of the anomalous dimensions 1)4.28)1 is different compared to the ones 
found in the previous section, Eq. ()4.8jl . The agreement is achieved by means of the Landen 
transformation of the modular parameters jTHj 



(i-vT 



-(0) 

^' = - V V!-^^ ■ (4 - 33) 

upon which the relations 1)4.28)1 and ()4.8jl coincide provided that m = 1. Remember that (— x) 
depends on the radial coordinate of the spike po, Eq. ()4.7jl . so that vl — r = e~ 2po and the long 
string limit po - >■ 00 corresponds to r — > 1. We have demonstrated in Sect. 12.31 that on the 
gauge theory side the limit r — > 1 corresponds to a — > l/(2m) and 6^0, Eq. (|3.13jl . As a 
consequence, the two cuts [—a, —b] and [b, a] collide at the origin yielding the logarithmic scaling 
of the one-loop anomalous dimension (J2.42)) . Let us examine the limit r — >• 1 of the obtained 
stringy expressions ()4.27)1 and ()4.28|) . 

Since K(r) ~ — |ln(l — r) for r — > 1, one would expect from ()4.27j) that 6 str should vanish 
in this limit. Indeed, this is the case for A' = while for A' 7^ one deduces from ()4.27jl that 
the reality condition for 6 str implies that its minimal value is bounded as 6 str > V^V. Carefully 
examining 1)4.27)1 for r — > 1 one finds 



-'str 



m 2 In 2 Vl — t ^ 



For A' = this relation coincides with the one-loop expression b^l = 2b/ (3, Eq. (|3.13|) . Matching 
(14.34)1 into ()4.29)1 we conclude that higher order corrections to 6 str push its minimal possible value 
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away from the origin and, therefore, prevent the two cuts [— a str , — b str ] and [6 s tr,Ostr] to collide. 
From ()4.28j) one finds in the limit r — > 1 

P-^r T = + • • • , 7(A) = L [-mb stx In - l] + . . . . (4.35) 

The asymptotic behavior of these expressions depends on the value of the parameter £ str = 
A'ln 2 (iV/L), Eq. (gUJ). 

For £ str < 1 the expression for 6 str , Eq. (|4.34|) . admits a series expansion in A' and leads 
together with (|4.35J) to 

1 N , LIN 

»str = : N + • • • , flstr = Tr—F + • • • ' Vl - r = — — ^ + . . . . 4.36 
m In — 2mL 2 In — 

Substitution of ()4.34|) into ()4.35|) yields an expression for the anomalous dimension 7(A) which 
coincides with ()4.20|) for m = 1. It is instructive to compare the positions of the cuts in gauge 
theory, a and b, and on the string side, a = f3a stI /2 and b = (3b stv /2. Here the additional factor 
(3/2 appears due to a different definition of the local parameter x on the spectral curves ()2.37|) 
and (JH22D- In this way we find a = I /(2m) + . . . and b = (L/N)/(4mln(N/L)) + ... which 
coincides with the similar relation (j3.13J) in gauge theory to one-loop order. 

For £ str ^> 1 the expression for b stT , Eq. (j4.34J) . is not analytical in the BMN coupling A' = 
A/(vrL) 2 

m Lin 
b str = VA 7 + . . . , a str = — — + . . . , Vl -r = — VA 7 + . . . , (4.37) 

2mL 2Jy 

while the expression for a str is the same as in ([4.36)1 . This suggests that for ^ str ^> 1 higher order 
corrections only modify the lower edge of the cut. One finds from ()4.35|) that the anomalous 
dimension scales as 7(A) ~ m\^X\n(N/(m\/X)) and matches (j4.21|) for m = 1. We conclude that 
the logarithmic scaling of the anomalous dimension for iV ^> L on the string side is realized when 
the inner boundary of the cut 6 str approaches its minimal possible value x/A 7 (see Fig. Efc) which 
coincides with the position of the double pole of the differential ()4.25[) . 



5 Conclusion 

In the present paper we have studied the properties of anomalous dimensions of Wilson operators 
of higher twist L and large Lorentz spin iV in the weak and strong coupling regimes by making use 
of the remarkable integrability symmetry on both sides of the gauge/string correspondence. We 
concentrated on operators which have the minimal anomalous dimension for given Lorentz spin 
and put a special emphasis on the appearance of the single-logarithmic behavior, 7(A) ~ IniV, 
in the thermodynamic limit L — > 00. 

On the gauge theory side, we applied the method of the Baxter Q— operator to identify 
different regimes of the minimal anomalous dimension in integrable sectors of (supersymmetric) 
Yang-Mills theory to one-loop order. We argued that for iV ^> L the asymptotic behavior of 7(A) 
is controlled by the parameter £ = ln(N/L)/L. For £ < 1 the anomalous dimension possesses 
the BMN scaling 7(A) ~ A/L, while for £ ^> 1 it scales logarithmically with N. Transition to the 
second, logarithmic regime manifests itself through the divergence of the semiclassical expansion 
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for 7(A) as £ ~ 1. The anomalous dimension is uniquely determined by the configuration of Bethe 
roots which condense in the thermodynamic limit on two symmetric cuts. We demonstrated that 
the semiclassical approach breaks down for £ ~ 1 due to the collision of cuts at the origin and 
worked out an asymptotic expression for anomalous dimensions which is valid throughout the 
entire region of £. 

On the string theory side, we used the identification of the minimal anomalous dimension of 
scalar operators in the Af = 4 SYM at strong coupling as the energy of folded string rotating on 
AdSsxS 1 part of the target space. Similar to the previous case, the anomalous dimension has 
different behavior for N ^> L depending on the value of the parameter £ str = X\n 2 (N/L)/L 2 . For 
£ str < 1 the anomalous dimension has a regular expansion in powers of the BMN coupling and 
its lowest term matches the one-loop expression for 7(A) at weak coupling for £ < 1. For £ str S> 1 
the anomalous dimension scales logarithmically but its dependence on the 't Hooft coupling is 
not analytical anymore. We described the latter regime using two different (although equivalent) 
configurations. In terms of the folded rotating string the logarithmic scaling occurs when two 
most distant points of the string (two spikes) approach the boundary of the AdS space. In terms 
of the spectral curve for the classical string sigma model, the same configuration is described by 
the elliptic curve with symmetric branching points. Different regimes of 7(A) arise depending 
on the position of the branching points. In the logarithmic regime, the inner branching points 
approach the minimal possible value ±vA' so that the anomalous dimension ceases to obey the 
BMN scaling. 

Integrability played a key role in our analysis. In generic (supersymmetric) Yang-Mills theory 
it holds to one-loop order for Wilson operators belonging to special, holomorphic sectors only. We 
would like to stress that logarithmic behavior of anomalous dimensions is not tied to integrability. 
In non-integrable sectors the mixing matrix for Wilson operators contains additional terms which 
break integrability symmetry. They do not affect however the logarithmic scaling of anomalous 
dimension for iV — > 00. The reason for this is that the logarithmic scaling of anomalous dimension 
can be associated with the contribution of soft gluons (i.e., gauge field quanta). Soft gluon 
radiation is not sensitive to the quantum numbers of Wilson operators (except of the total 
Lorentz spin) and, therefore, it provides the same logarithmic contribution to the anomalous 
dimension of Wilson operators in all sectors. Integrability allows one to identify various regimes 
of the asymptotic behavior of anomalous dimensions and to determine the "critical" value of 
\n(N/L) ~ L at which the logarithmic scaling sets in. 

As a function of the total Lorentz spin N, the anomalous dimensions of twist— L operators 
occupy a band. Our discussion was restricted to the minimal anomalous dimensions belonging to 
the lower edge of the band. It would be interesting to extend the above analysis to excited states 
(13.47(1 with m > 1 and describe the band structure occupied by the anomalous dimensions as it 
arises form the Baxter equation on the gauge theory side. In string theory, there are two different 
classical configurations yielding the same logarithmic asymptotics of anomalous dimensions at 
strong coupling — the multiple folded string ^T] and the spiky string [38 j. The coefficient in 
front of the logarithm, m, is twice the number of foldings of the string on itself, in the former 
case and it is the number of spikes, in the latter one. It would be interesting to construct a 
generic configuration which interpolates between both solutions. 

Note added: After the work has been completed we were informed by Yuji Satoh that he, 
in collaboration with Kazuhiro Sakai, analyzed the spectrum of anomalous dimensions with a 
special emphasis on the appearance of logarithmic scaling in the large spin limit. For the two- 
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cut solution, their findings are in agreement with the analysis presented in sections 13.11 and 14.21 
of this paper. In particular, they also came to the conclusion that, in the thermodynamical 
limit L — > oo, the semi classical approach to the Bethe Ansatz equations is only applicable for 
]n(N/L) < L and one can not reproduce ~ IniV behavior of anomalous dimensions unless the 
finite size corrections in 1/L are included. 
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Appendix Asymptotic expression for the energy 



The asymptotic expression for the energy ()3.37j) involves roots of the transfer matrix and is not 
well suited for performing the thermodynamical limit L —>■ oo. Let us rewrite the energy in terms 
of the transfer matrix itself. We notice from (j3.30J) that 



d_ 

du 



L 1 
< • ii — 



k=l 



u — 5k 



(a.i; 



and replace the ip— functions in ()3.37j) by series representation 

oo 

tp(s + i5 k ) + ip(s - i5 k ) - 2ip(2s) = £ 



k=0 



+ 



— I 



i(s + k) — 5k —i(s + k) — 5 k 2s + k 



(A.2) 



Then, the energy ()3.37j) can be written as an infinite sum 



d 



£7 = 21n2 -£)—/(*), /(x) = ln 



k=0 



t L (i(s + x))t L (-i(s + x)) 



(2s + x) 



2L 



(A.3) 



where f{x) — > 2 In 2 for x — ► oo. The sum can be evaluated with a help of the Euler-Maclaurin 
summation formula 



1 _ R A 

e = /(o) - \f'{0) + £ pyy/ (2fc) (0) = 



e d -l' 



(A.4) 



with B 2 k being Bernoulli numbers. Taking u = ±i(s + x) in the Baxter equation ()2.3|) . one 
obtains ti(±i(s + x)) as the ratio of Q— functions and finds for x — > 



f(x) = In 



Q(i(s + x + 1)) Q(-i(s + x + 1)) 
Q(i(s + x)) Q(-i(s + x)) 



0(x L ) . 



(A.5) 
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Then, one uses the WKB ansatz Q(u) ~ exp(r/ 1 S(ur])) to get 
f(x) = ^[S(i(s + x + l)r])-S(i(s + x)r))+S(-i(s + x + l)r))-S(-i(s + x)r])]+0(x L ) . (A.6) 

Substituting this expression into (jA.4j) one finally obtains the relation 

1 26 
E = /(0) - -/' (0) + . . . = -2s V S"(0) + ... = -fS"(0) + . . . , (A.7) 

which coincides with the semiclassical expression ()2.33|) . It is important to keep in mind that 
Eq. (jA.7|) was obtained under the assumption that contribution of terms with higher derivatives 
is small, /"(0) <C /'(0), or equivalently r}S"'(Q) <C S"(0). This relation is not satisfied if the 
distribution of Bethe roots scales at the origin as ~ Inx, Eq. (|2.4U|) . In that case, one has to rely 
on the formula (jA.3|) which resums all singular terms. 
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